NASA TECHNICAL. NOTE 



NASA TN D-7845 


LOAN COPY: RETw 
AFWL TECHNICAL £ 
KIRTLAND AFB, Rj 


CALCULATION 




WASHINGTON, D. C. 


TECH LIBRARY KAFB, NM 



TECH LIBRARY KAFB, NM 



□ 133b22 


1. Report No. I 2. Government Accession No. 

NASA TN D-7845 j 

4. Title and Subtitle 

A COMPUTER PROGRAM FOR CALCULATION OF SPECTRAL 
RADIATIVE PROPERTIES OF GAS MIXTURES 

7. Author(s) 

John E. Nealy 


3. Recipient's Catalog No. 


5. Report Date 

August 1975 


6. Performing Organization Code 


8. Performing Organization Report No. 

L-9745 


9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, Va. 23665 


10. Work Unit No. 

502-27-01-03 


11. Contract or Grant No. 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

15. Supplementary Notes 


16. Abstract 


13. Type of Report and Period Covered 

Technical Note 


14. Sponsoring Agency Code 


A computer code is described whereby calculations of radiative properties of gas 
mixtures may be made. The program is arranged so that distinct radiative processes for 
each species are computed in individual subroutines. Provision is made for calculating 
radiative properties in nonequilibrium situations - separate rotational, vibrational, and 
electronic temperatures may be used. These features should provide a flexibility not 
currently available in such programs. 

The basic equations and the program listing in FORTRAN IV language are presented. 
Sample calculations are given for high -temperature air and carbon dioxide and are compared 
to calculations made with previously developed programs. 


17. Key Words (Suggested by Author (s) ) 


18. Distribution Statement 


Thermodynamics 
Spectroscopy 
Thermal radiation 
Radiative heat transfer 


Unclassified - Unlimited 


Subject Category 33 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* 

Unclassified 

Unclassified 

86 

$4.75 


For sale by the National Technical Information Service, Springfield, Virginia 22161 



A COMPUTER PROGRAM FOR CALCULATION OF SPECTRAL 
RADIATIVE PROPERTIES OF GAS MIXTURES 

John E. Nealy 
Langley Research Center 

SUMMARY 

A computer code is described whereby calculations of radiative properties of gas 
mixtures may be made. The program is arranged so that distinct radiative processes 
for each species are computed in individual subroutines. Provision is made for calcu- 
lating radiative properties in nonequilibrium situations - separate rotational, vibra- 
tional, and electronic temperatures may be used. These features should provide a 
flexibility not currently available in such programs. 

The basic equations and the program listing in FORTRAN IV language are pre- 
sented. Sample calculations are given for high -temperature air and carbon dioxide and 
are compared to calculations made with previously developed programs. 

INTRODUCTION 

Radiative properties of gases are of interest for a variety of reasons, and a great 
number of complex and laborious experiments and calculations have been made on the 
subject. Basic radiative processes of gas constituents are of prime importance for 
application to energy -transfer problems such as occur in stellar and planetary atmos- 
pheres, combustion reactions, shock layers of reentry vehicles, and high-temperature 
plasmas. 

The program described in this paper has been developed because of a continuing 
need for a more versatile and general code which can be used to predict the spectral 
radiative properties of a mixture of gases. Some programs for calculating such prop- 
erties are in existence (refs. 1 to 5) and have been usefully applied to a number of 
radiative transfer problems (refs. 6 and 7). For example, the program of reference 1 
has been written primarily for application to radiative flow-field problems at the sacri- 
fice of spectral detail. The program of reference 2, on the other hand, provides con- 
siderable spectral detail but, because of the large storage requirements, has had more 
limited application to radiating flow fields. In addition, it was originally developed for 
radiating species of air, and generalization to arbitrary mixtures is somewhat involved. 



This program has been developed in part from earlier programs written by the 
author for calculating specific radiative processes. A separate smearea-rotational- 
line band system program has been exercised previously on several occasions (ref. 8) 
and is included here as a subroutine. Other distinct processes are also handled by 
individual subroutines. It is felt that such an arrangement adds greatly to the versa- 
tility of the program in that it facilitates inclusion of modifications to the various com- 
putational procedures. 

It is apparent from the literature dealing with gas radiative properties that basic 
input parameters require updating at frequent intervals. Such input changes should be 
straightforward and easily accomplished if a radiative-properties program is to retain 
its applicability. Throughout the development of the program described herein, consid- 
erable effort has been directed toward keeping the input simple and well defined. Since 
the object of the program is to compute radiative properties for a desired gas mixture, 
the radiative input data for an individual species comprise a completely independent 
unit containing all parameters required to calculate the radiance for that particular 
species. This feature greatly facilitates the inclusion or deletion of desired species 
from case to case and permits the existence of a readily usable punched-card input data 
file. 
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SYMBOLS 

Einstein spontaneous transition probability, sec-1 

rotational spectroscopic constant, cm - l 

Planck blackbody radiation function, W/m^-Mm-sr 

speed of light, m/sec 

energy level value, cm"l 

electronic charge, C 

oscillator strength, or f-number 

degeneracy 

Planck constant, J-sec 
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I ionization energy, cm~l 

Q 

j radiative volume emission coefficient, W/m -pm-sr 

k Boltzmann constant, J/K 

L path length, m 

M particle molecular weight, kg/kg- mole 

m electron rest mass, kg 

N particle number density, m _ 3 

Nq Avogadro number, m~3 

n principal quantum number 

P pressure, N/m^ 

Q partition function 

q Franck- Condon factor 

p 

R radiant intensity, W/m -pm-sr 

Ry Rydberg constant 

s parameter in empirical formula (19) 

T temperature, K 

v* upper vibrational quantum number 

v" lower vibrational quantum number 

W integrated spectral radiant intensity, W/m^-sr 

w line half -width, m 
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z 


degree of ionization 


a polarizability, cm3 

a 1 spectroscopic constant for vibration- rotation interaction, cm-1 

/3 parameter in empirical function (19) 

y partition function parameter defined by equation (42), cm - * 

6q permittivity of free space 

X wavelength, m 

ju absorption coefficient, m“l 

v frequency, sec"* 

v wave number, m"l 

cr cross section, m2 

p density, kg/m3 

P 0 density of air at standard temperature and pressure, kg/m 3 

w e vibrational spectroscopic constant, cm - * 

&> e x e vibrational spectroscopic constant, first anharmonic correction, cm-1 

cu e y e vibrational spectroscopic constant, second anharmonic correction, cm-1 

Subscripts: 

A general particle 

D Doppler broadening 

e electronic 
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eff effective value 

FB free-bound process 

FF free-free process 

H hydrogen atom 

L individual atomic line 

S. lower energy level 

N general neutral particle 

R resonance broadening 

r rotational 

S Stark broadening 

th threshold value 

u upper energy level 

v vibrational 

v' upper vibrational quantum number 

v" lower vibrational quantum number 

0 initial value 

COMPUTATIONAL OUTLINE AND GENERAL 
PROGRAM DESCRIPTION 

Basic radiative properties of various substances are well defined in the literature 
sources treating the subject. There is, however, wide diversity and nonuniversality 
of nomenclature existing in this field. The important basic properties used in this 
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paper will all be defined, therefore, and the governing equations used in the calculations 
will be developed with care. 

The spectral intensity of the ultimate radiant emitter-absorber is described by 
the Planck blackbody function for specific intensity 


bJt) = 


2hc2 

X5 


1 

e hcAkT _ t 


( 1 ) 


The linear absorption coefficient p is defined by the Lambert law of exponential attenua- 
tion of incident radiation Rq over path length L 

R = R 0 e _/iL (2) 

For an extended medium of absorption coefficient p , which also emits radiation with 
volume emission coefficient j, Kirchhoff's law may be expressed as 

j = PB°(T) (3) 

for media in local thermodynamic equilibrium. At a given condition, this equation may 
be written for a single radiative process; that is, j and p for various processes are 
additive quantities at a given wavelength. 

Another basic radiative property is the cross section for absorption, which is 
simply the absorption coefficient divided by the appropriate contributing particle-number 
density: 


"a'Na 


(4) 


In this program, p and a for the desired processes are computed as functions of 
wavelength and are used to obtain a value of j. The j's are then summed for all 
species included to arrive at the overall spectral emission coefficient for the mixture. 


For application to practical problems, the volume emission coefficient of an 
extended medium must be used in a radiative transfer equation which is integrated over 
the geometrical dimensions of the medium. The radiative transfer equation is generally 
difficult to solve, and only one specific case is dealt with here. This is the case of the 
one-dimensional slab of width L which extends to arbitrarily large dimensions in 
the other orthogonal directions. It may be shown that the emergent radiant intensity 
normal to the surface is given by 
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(5) 


R = B°(T)(l - e ^ 

For an optically thin medium, that is, for /i L« 1, 

^ = pB? = j (6) 

L x J 

To obtain the total energy emergent from a surface of the slab equation (5) is integrated 
over all wavelengths: 

W = J B°(T)(l - e _fxL )dX (7) 


This equation has been evaluated for several cases and compared to similar calculations 
performed with other programs. 

The radiative processes for each species (free-free, free-bound, lines, etc.) are 
also handled separately in the program. Various models for computation of these 
processes exist, and those chosen for use in this program are described in detail in the 
text. Empirical formulas may often be employed advantageously because of the compu- 
tational simplification achieved. Wherever these have been used in the program, the 
procedure is clearly outlined by which the empirical formulas are obtained from more 
detailed quantum-mechanical calculations. 

Radiative transfer problems often occur in situations for which the various con- 
tributing processes take place at different effective temperatures. Examples of such 
situations are the zone of chemical nonequilibrium in shock layers and plasmas for 
which the ion and electron temperatures are different. For each of the processes 
included in the program, an option exists of having a different electron temperature as 
input. For the molecular band systems, electronic, vibrational, and rotational temper- 
atures may be employed. 

An indication of the computational procedure is shown in figure 1, where an abbrevi- 
ated flow chart is presented. It is seen that the program consists of three sections: 

(1) The main program which handles input data and calls the various subroutines. 

(2) The group of subroutines which computes spectral radiative properties for each 
process. 

(3) The output section which accumulates and stores the various spectral properties 
as they are being computed and specifies the form of the final output. 
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Intermediate output for calculations of the specific radiative processes may also be 
obtained. 

A program listing and detail on the FORTRAN coding is given after a description of 
the equations used for calculating the various radiative processes occurring in the follow- 
ing subroutines : 

FFCON - computes free-free continua for all species 

FBCON - computes free -bound continua for atoms and atomic ions 

LINES - computes line radiation for atoms and atomic ions 

SMEAR - computes molecular band radiation for molecules with smeared-rotational- 
line model 

INFRA - computes infrared vibration-rotation spectra for molecules 
RACOM - accumulates, sorts, and adds spectral emission and wavelength points 
PUTTY - output section of program 

These processes are included in the sample calculations given later. Subroutines for 
other processes may be easily incorporated once a desirable computational procedure 
has been determined. 

CALCULATION OF FREE- FREE RADIATION 


Free-free radiation arises as a result of the process 
A+e - — A+e' + hi' 


in which a free electron e" interacts with the coulomb field of another particle A. 

The consequent change in energy of the electron results in an emitted photon. The 
cross section for the process is normally evaluated in units of (length)^ which, when 
multiplied by the number density product N A N g , yields the linear absorption coefficient 
in reciprocal length units. Detailed quantum-mechanical calculations of these cross 
sections may be reliably made for relatively simple atoms and molecules. An excellent 
procedural outline on such calculations is given in the text by Armstrong and Nicholls 
(ref. 9). 

The classical expression for the emission cross section of a free electron accel- 
erated in the field of a particle of charge Z was developed by Kramers and may be 
written (see, for example, ref. 10): 


16t 7 2 e 6 A 3 Z 2 1 

hc4(2>rm e ) 3/2 < kT > 1/2 (im 0 f 


( 8 ) 
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The corresponding emission coefficient resulting from the use of equations (1) and (3) is 


. 32,2 Z 2 e 6 N e N A e*p(-J^) 

3 l/3 (47re 0 ) 3 c 2 (2 7 rm e ) 3 / 2 (kT) 1 / 2 X 2 


with the factor 1 - exp(-hc/xkT) included to account for induced emission. For elec- 
trons interacting with ions, this expression appears to give reasonable results (see 
ref. 11), and has been used in the program in this form. 

For electrons interacting with neutral particles, a first approach is to measure or 
deduce an effective Z 2 for the particle in question and to apply the Kramers equation 
(eq. (9)). This procedure may be used to give rough order -of -magnitude cross sections 
for a limited temperature and frequency range. In reference 3 such effective Z 2 
parameters have been used over an extensive wavelength and temperature range, but the 
resultant inaccuracies are minor since free-free radiation from neutral particles is not 
a significant contributing process in the overall radiance of high-temperature equilibrium 
gas mixtures. 

Detailed calculations of neutral Bremsstrahlung (refs. 12 and 13) indicate that an 
effective Z 2 must be a function of X and T. The cross-section calculations of 
reference 12 have been used here to evaluate an approximate wavelength and temperature 
dependence of such an effective Z 2 . This variation has then been incorporated into 
equation (9) with Z 2 * Z 2 ff , which has been used in the neutral free-free radiation 
contributions of the program. The term Z 2 ££ is determined from the empirical 
formula which fits the calculations of references 12 and 13: 


r7 2 _ n2 3(10)~ 10 (XT - 20 000) 2 

Z eff " Z N -^572 


+ t 5 / 2 (10)' 10 


(10) 


where X is in units of pm and Z?. is now an empirical constant. Since equa- 

^ 2 2 

tion (10) only approximates the X,T dependence of Z e ff, Zjq- has been evaluated at 
T = 10 000 K, X = 2 pm in order to match cross-section data of references 12 and 13. 


Where detailed cross-section calculations are not available, another approach has 
been utilized. An approximate free-free absorption cross-section formula has been 
derived in reference 14 for highly polarizable neutral atoms. The polarizability a 
of a neutral particle has been used to estimate the free-free absorption cross section as 


a FF = 1.62(10) -31 \fcx 3 T 


(ID 
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Comparison with equation (8) indicates that Z^ff should be given by 

Z^ ff = 11 900T 3 / 2 /c7 (12) 

(The constants in eqs. (11) and (12) apply when X is in pm and a is in cm3.) The 
accuracy of equation (11) is best at high temperatures (above about 10 000 K). Its appli- 
cability also decreases as a decreases. However, the values of z|j£ determined 
from equation (10) compare rather favorably with those deduced from equation (12) 

(at T = 10 000 K), and use has been made of equation (12) when detailed cross-section 

2 

calculations have not been available. Table I gives the Zjq- used for all neutral 
particles in the illustrative calculations presented later. A comparison of ov,t- 
determined from equations (8) and (10) with the cross sections of reference 12 for the 
carbon atom is given in table II. The agreement over the range of T and X is 
typical for the other species considered thus far. 

CALCULATION OF FREE-BOUND RADIATION 

As its name implies, free-bound radiation occurs when a free electron combines 
with an ion and excess energy is emitted as a photon; thus 

A + + e — A + hv 


The absorption cross section (for the inverse process) may be expressed by a 
relatively simple equation for the hydrogen atom (see, for example, ref. 10): 


ct FB,H 


64t r 4 
3/3 


e 10 m Z 4 X 3 
e 

h®c 4 n 3 ^47reQ^ 


(13) 


In order to account for the nonhydrogenic species, several factors should be taken 
into consideration. For the absorption process, only allowed transitions between the 
parent atom and its ion should be included. In addition, the hydrogenic formula (13) 
is inaccurate for transitions involving the lower lying states of the parent atom. It 
becomes increasingly more applicable for higher levels, or increasing values of n. In 
the present version of the program equation (13) is used for n ^ 3. 

O 

As is done in reference 15, an effective Z is defined by 


7 2 

z eff,FB = 


2 / E u - E f 

\ R Y 


(14) 
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The cross section then becomes 


647r 4 el0m e Z 4 A 3 
cr(n) = £ 

3\/3(47re 0 ) 5 h 6 c4n 5 

For the hydrogen-like particle of ionization energy I, 



(15) 


647r 4 e 10 m e A 3 I 2 
11 3v / 3(4 7 re 0 ) 5 h 6 c 4 R|.n 5 

The absorption coefficient is found by summing the product of a(n) and the particle 
number density in level n: 


V = 


2N1 2 

Q 


647r 4 e 10 m e A 3 

3\[3(47rA 0 ) 5 h6c 4 R|. 


I 


4jexp 

n' 3 


I 

n 2kT 


(17) 


_ g 

M = / N(n) cr(n) = — / g(n) a(n) exp — - (18) 

U Q L-i kT 

n n 

where the approximate hydrogenic upper-level degeneracy, g(n) = 2n^ has been used. 

As with the free-free processes, detailed quantum-theory calculations for a 
variety of simple species are available, and for computational convenience some measure 
of empiricism is desirable. It is particularly important to account for nonclassical 
effects for the lower lying levels of nonhydrogenic species. The deviation from hydro- 
genic cross sections is most severe for these levels, and for high-temperature transfer 
calculations transitions involving these levels are among the most important contrib- 
utors. A relatively large number of such free -bound cross sections are reported in the 
literature and have been calculated in various degrees of detail. Due to the importance 
of N and O atoms in atmospheric photochemical processes and shock layers of reentry 
vehicles, their free -bound behavior has been studied rather extensively, both experimen- 
tally and theoretically (see refs. 16, 17, and 18). More recently, the determination that 
Mars and Venus have atmospheres composed primarily of CO 2 has led to increased study 
of the radiative properties of carbon (see refs. 18, 19, and 20). 
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A large number of ground-state connected free-bound cross sections are reported 
in reference 21. Each of the cross sections is represented by an empirical formula: 


ct FB " 



X£X th 


(19) 


The parameters s and j3 have been determined so that a close fit between equation (19) 
and a detailed wave-function calculation is obtained. The function (19) is used in the pres- 
ent version of the program, and values of parameters s and /3 as given in reference 21 
for C, N, and O atoms are used as input. Calculations of apg for these atoms using 
equations (15) and (19) in combination are given in table III. 

CALCULATION OF LINE RADIATION 


The radiation emitted when electronically excited atoms undergo transitions to 
lower levels gives rise to line radiation of wavelength X^ u , in which the photon energy 
of such a transition is 

^-*E u -E £ (20) 

x Hu 


The radiant intensity of the line is given by 


!u( = 5; N >A.t ^ 


( 21 ) 


where A uJ ^ is the Einstein coefficient or spontaneous transition probability. The rela- 
tionship between the spontaneous transition probability and the absorption oscillator 
strength, or f-number, is (see, for example, ref. 22): 


Since 


\lS. = 


H 87r 2 e 2 f lu 
£u m g cA^ 47^0 


N f S £ X fu kT 


( 22 ) 


(23) 
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and 


N 


-Eje 


N * = ex P^ 


kT 


(24) 


the volume emission coefficient for an individual atomic line may be written as 




2he' : 


( 47 re o) m e 


N gjg f fu 
Q .3 
a £u 


exp 


-< E i * 


he) 


4u 


kT 


(25) 


This equation is used in the present calculations for line radiancies. The expression 
gives only the total power radiated in the transition, with no information about the 
spectral profile of the line. 

Spectral-line profiles are almost entirely dependent on the state of the medium in 
which the line is emitted (p, T, N e , etc.). The inclusion of the various line -broadening 
mechanisms into a prediction of half-widths for polyelectronic atoms is very complicated, 
usually accompanied with a significant degree of uncertainty. In high-temperature 
plasmas, however, the energy contained in the wings of the lines is of importance and 
broadening must be accounted for. Several good state-of-the-art approaches to this 
problem are described in a number of sources. See, for example, references 9, 10, 22, 
and 23. These references indicate that a function describing combined Doppler and 
collision broadening is best for the majority of applications. One such function is 
commonly referred to as the Voigt profile (ref. 9) and is relatively cumbersome math- 
ematically. For this reason, it is often used for line-shape calculations in one of its 
asymptotic forms. Even though the Voigt function has not been incorporated into this 
program, it should be recognized that circumstances occur which warrant its use. 

The collision, or Lorentz, profile is used for line calculations presented herein. 

It is particularly applicable to Stark broadening, which dominates in low-density ionized 
gases. As the degree of ionization in a gas decreases, the Lorentz profile becomes more 
conservative; that is, there is generally an overprediction of the radiant energy con- 
tained in the line wings. 

The spectral emission coefficient for a Lorentz line is given by 



W L 

(A - A UJ? ) 2 + w^ 


Using Kirchhoff ’s law, the absorption coefficient may be written as 


(26) 
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(27) 


M ^(x) 


3L W L 

iB (V)( x - x uj ) 2 + w l 


The half -width of the line for Stark broadening w Lg is determined in this paper 
by the approximate formula described in reference 24: 


w Ls = 5 - 7 ( 10 ) 


-22- N e 4f f T- l/2 

X 2 (Z + l) 2 


(28) 


where the proportionality constant is valid for SI units. When the effective quantum 
number of the lower bound state n e ff is defined as 



(Z + l) 4 R 2 

(I - E *) 2 


(29) 


the line half -width formula becomes 


w 


LS 


5.7(10)" 22 


(Z + l) 2 X 2 N e R^ 

T 1 / 2 )! - E £ ) 2 


(30) 


Doppler broadening may contribute significantly at high temperatures and densities. 
Pure Doppler broadening results in a line shape with wings decaying exponentially instead 
of quadratically, as is found for the Lorentz shape. The half -width for Doppler broadening 
is given by 


W LD=*' 


2 N Q kT log e 2 
Me 2 


(31) 


Resonance broadening is important for strong lines of a species in a mixture where 
this species has a large number density. Collisions between these identical particles 
have a much larger effect on line half -width than do collisions between unlike particles. 
The relationship for line half-width for this mechanism is (ref. 23) 


W LR = 


3tt g£ N ffu e2x3 
^47re 0 ) gu 27rm e c^ 


(32) 


14 



I 


where N is the number density of the particular species involved. The line shape for 
resonance broadening is Lorentzian. 

Equations (25) and (26) have been used to compute line radiation from atoms and 
ions in all illustrative calculations presented later. The line half-width has been taken 
to be the sum of the half-widths of the Stark and resonance broadening processes pre- 
viously described 

W L = W LS + W LR + W LD (33) 

Where Doppler broadening is a significant mechanism, however, equation (33) applied 
to a Lorentzian profile would be inaccurate. In such circumstances, a Doppler or Voigt 
profile would be more appropriate. For the calculations presented here, and for a wide 
range of pressure and temperature, equation (33) provides sufficient accuracy. 

All basic input data for atomic and ionic lines have been taken from reference 25. 
Included in the sample calculations are line contributions from the species C, C + , C ++ , 
N, N + , O, and 0 + . 


CALCULATION OF DIATOMIC MOLECULAR 
ELECTRONIC BAND SPECTRA 


The radiation resulting from transitions between electronic states in molecules is 
much more complex than the corresponding atomic process. This complexity is due to 
the additional degrees of freedom by rotation and vibration of the atoms composing the 
molecule. Instead of a single line, a many-line spectrum or band system results and 
is seen in a generally well-defined pattern over an extensive wavelength region. The 
position of the individual lines of a band system is determined by the vibrational and 
rotational levels within the two electronic states involved in the transition. The spacing 
between individual lines arising from adjacent rotational levels is usually very narrow, 
and a commonly used spectral model has been developed which treats these closely spaced 
lines as a continuum. This model is usually referred to in the literature as the "smeared- 
line" model, and detailed derivations may be found in references 26 and 27. 

The positions of the various band heads of a band system are determined by the 
vibrational levels of the two electronic states 


> - he 

V ’ V "~ 


(34) 
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where 


E u (v ') = E u (0) + he 


2 , 
c ° e ( V ’ + I) " c °e x e( v ’ + l) +w ey e ( v ’ + |) 


(35) 


and 


E £ (v") = E £ (0) + he 


“e( v ” + |) - “e x e( v " + if + “ e ye( T " + if 


(36) 


The spectral absorption coefficient for a molecular band system as derived in 
references 26 and 27 and modified to include separate electronic, vibrational, and rota- 
tional temperatures as has been done in reference 28 is given by 


MOO = ^ N, ^0 f 


B, 




me 2 * * ^ u Q v k £ T r 7re 0 


(37) 


Here all energy terms and spectroscopic constants are considered to be in units of 
reciprocal length. The term Xqq refers to the band-head wavelength for v’ = v" = 0. 
The function <p gives the spectral variation throughout the band system resulting from 
the v" - v* transitions 


-II 


^v , v t 


v' v" l B v' 


- B, 


exp 


E £ (v") 

kT„ 


hcB v” /l . 1 \ 

k \X X V ' v"/ 

T (B , - B 
V\ V' v'7 


(38) 


The summation extends over the range of upper and lower vibrational quantum numbers 
and is nonzero only for (a) X >X v , v „ when the band is red-shaded, and (b) X <X v , v „ 
when the band is violet -shaded. The Franck-Condon factors q y t v n are evaluated from 
vibrational-state wave functions and a large number of sources exist in which these 
arrays are tabulated. For the present calculations, values listed in reference 28 have 
been used when possible. The quantities By appearing in equation (38) are found from 





The electronic absorption oscillator strength for the band system f^ u is defined 
in reference 10 and is often assumed to be wavelength independent. As is pointed out in 
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references 29 and 30, such is not the case, and provision for a simple functional variation 
has been anticipated. Since the calculations appearing here have been carried out for 
purposes of comparison, however, no X-dependence of f^ u has been employed. 

In order to evaluate the number density of particles in the absorbing state 
the partition function formula of reference 31 has been used 




y Q(n) 


n 


where 




exp 


hew 

"kT, 


v-1 


eSL 


and 


n 


1 too „ a a \u 

_J_ + J_ + 8 _L_ he 

w e,£\ w e,£ B f w e k 


(40) 


(41) 


(42) 


In the preceding equations, the form of the partition function has been slightly modified 
in order to account for different electronic, vibrational, and rotational temperatures. 

The diatomic electronic band systems included in the present calculations are 
given in table IV with references to the input data for each band system. The data are 
presented later in a punched-card listing. Many band-system calculations have been 
made with this subroutine and appear to agree well with similar computations reported 
in references 26, 28, and 32. An example of the subroutine output plotted for the CO + 
comet tail band system is given in figure 2 where the effects of different rotational and 
vibrational temperatures are illustrated. 

CALCULATION OF RADIANCE FROM INFRARED VIBRATION-ROTATION 
BANDS OF DIATOMIC MOLECULES 

At moderate and low temperatures (below about 5000 K) the infrared emission- 
absorption properties of gases assume greater importance. In an emitting gas, the 
fractional portion of the radiant energy contained in the infrared generally increases with 
decreasing temperature, and molecular gases are less dissociated at lower temperatures. 
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Infrared vibration- rotation spectra arise as a result of transitions between vibra- 
tional levels within a given electronic state. This electronic state is commonly the ground 
state, since the pertinent temperatures are usually so low that a great majority of the 
molecules are in this state. In the calculations of this paper, only two diatomic species 
which are not homopolar are of importance (NO and CO) and results are presented for 
these. (Homopolar molecules such as N 2 , 0 2 , H 2 , etc., have no vibrational spectra due 
to the absence of a dipole moment.) The equations used should apply equally well for 
other species if adequate input data are available. 

A model based on overlapping lines for these infrared bands has been developed by 
Breene (ref. 33) and is used in this program. The spectral absorption coefficient is 
given by 


^vV'W = 


_ -hc(A?) 2 

-E v * exp — — 

o N exp - — 4B f kT . . 

££_ kT v exp -he (A v) L~ f 

2 Q „„2 P 2kT 1 1 v'v" 
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= x - 

A Vv" 
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(44) 


Here f v , v , t is the oscillator strength for a given v* - v" transition and is 

the corresponding wavelength at the band center. The vibrational energy levels are 



(45) 


and the band-center wavelengths may be found from 
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(46) 


The quantities B y , are computed from equation (39), and the partition function Q is 
given by 


Q 



kT 

hcB 


(47) 


The individual band oscillator strengths may be computed from vibrational wave 
functions (refs. 34 and 35) or may be measured under conditions for which only one 
vibrational transition of the ground state is involved (ref. 10). The relationship between 
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ground-state vibrational transitions and those involving higher levels may then be used to 
determine other oscillator strengths. In reference 10, the ratios ^ v i v "/^iq are tabulated 
for the quantum- mechanical harmonic oscillator. For the present calculations, computed 
oscillator strengths for CO (ref. 35) and for NO (ref. 34) have been used. Equation (43) 
has been employed in the computation of the spectral profiles for the fundamental 
Av = 1 and the first and second overtones Av = 2,3 for these two species. 

In figure 3 one of the NO calculations of Breene (ref. 33) has been essentially 
repeated. Minor differences appear, probably because fewer total vibrational levels have 
been included, and slight differences in partition -function calculation may exist. In 
figure 4 the present-program computation for CO is compared to the more refined cal- 
culation of reference 36. The agreement is shown to be good in the region of high 
absorption. 


DESCRIPTION OF INPUT AND PROGRAM LISTING 

The principal function of the main program is to define the case to be computed, 
read the input parameters for each species, and call the subroutines related to each 
particular radiance contribution. How this is accomplished is most readily seen from 
the main program listing (appendix A) and the description of the input deck setup and 
key variables. A schematic of a typical deck setup is shown in figure 5. 

The first card to be read contains the following: 

Card 1: [NSPEC, TT, ELEC Format (13, 2E14.6)] 

NSPEC - total number of radiating species (does not include electrons) 

TT - particle translational temperature, K 
ELEC - electron number density, cm"3 

At this time a counter NSP is defined and is increased by 1 for each species encountered. 

The data input which follows is in the form of subdecks, each relating to a partic- 
ular species. These species decks are independent and may be arranged in any order, 
with the following exception: the first species should be an atom for which a free-free 
contribution is computed. The minimum of the wavelength range for this free -free must 
be the absolute minimum wavelength encountered in the remainder of the calculations. 

If it is not desired to include atomic species, input for a dummy species which does not 
contribute may be used to establish the minimum wavelength. Examples of dummy inputs 
may be found in the species input list (appendix B). 

The subdecks relating to atomic or ionic species contain input for the free-free, 
free-bound, and line contributions, respectively. The first card in a species deck defines 
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the species and specifies its number density: 

Card 1A [SPEC, ISPEC, NRAD, CON Format (A6, 213, E14.6)] 

SPEC - species identifier (chemical symbol) 

ISPEC - code distinguishing molecules from atoms (0 for atoms, 1 for 
molecules) 

NRAD - number of radiative contributions for this species (for example, if 
only free-free and free-bound are to be calculated, NRAD = 2) 

_ Q 

CON - species number density, cm 

For each contributing radiative process, a label card is read: 

Card FFl |TYRAD Format (7A10)] 

TYRAD - describes type of process being computed (for example, 

CARBON ATOM LINE RADIATION) 

For atomic or ionic species, the free-free input is read next: 

Card FF2 [NEQ, NW, WMIN, DELW, ZZ Format (12, 15, 3E14.6)] 

NEQ - specifies whether this process is to be computed for a temperature 
other than the local translational temperature (NEQ = 0 for equilib- 
rium calculations; otherwise its value is 1) 

NW - number of wavelength points at which the free-free process is 
computed 

DELW- increment between wavelength points, pm 
ZZ - appearing in equation (10) 

At this time subroutine FFCON is called and the radiance calculation is made. If 
NEQ = 1, an electron temperature, TTRON, Format (E14.6), is read and used in the 
calculation. The volume emission coefficient is computed for the specified number of 
wavelengths NW beginning at wavelength WMIN. In addition to the wavelength and 
emission coefficient, the cross section, absorption coefficient, blackbody function, and a 
running integral of the radiance are also included in the output. 

If this is the first calculation, the wavelengths are placed in an array WF(L) and 
the corresponding emission coefficients in F(L). This procedure is carried out in order 
to establish the radiance versus wavelength function to which the remaining contributions 
are added as the calculation proceeds. In the remainder of the calculations, the spectral 
radiance functions for a specific process are defined as G(L) and WG(L). For each 
subsequent radiative process, G(L) and WG(L) are added to the existing function F(L) and 
WF(L) to yield a newF, WF array. This addition of each individual contribution to the 
accumulated spectral distribution is performed in subroutine RACOM. 
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For an atomic species, the free-bound contribution is computed next. (FORTRAN 
variables common to the various subroutines which have been previously defined will not 
be redefined.) 

Card FBI [TYRAD Format (7A10)] 

Card FB2 (NEQ, NSTAT, NEDGE, XMIM, DELX 

Format (313, 2E14.6)] 

NSTAT - the number of states included in the calculation of partition function 

NEDGE - number of edges (thresholds for free-bound contributions) 

XMIM - minimum wavelength, pm (must be greater than WMIN in first 
FBCON input) 

DELX - wavelength increment establishing threshold discontinuity, pm 
(here always taken as 0.0001 pm) 

Cards FB3 . . . |GA(I), EA(I), I = 1, NSTAT Format (F4.0, E14.6)] 

GA(I) - degeneracy of state used in partition function calculation 

EA(I) - energy of state used in partition function (units of equivalent degrees 
K) 

Cards FB4 . . . |GL(I), PQN(I), EL(I), EU(I), I = 1, NEDGE 

Format (F4.0, 3E14.6)] 

GL(I) - degeneracy of initial state 

PQN(I) - principal quantum number of initial state 

EL(I) - energy of initial state (equivalent K) 

EU(I) - energy of final state, K (for the hydrogenic approximation, equa- 
tion (18), EU(I) is the ionization energy of the particle) 

Cards FB5 . . . [eSS(I), AL(I), SIG(I), I = 1, NEDGE 

Format (3F7.3)] 

ESS(I) - parameter s in equation (19) 

AL(I) - parameter £ in equation (19) 

SIG(I) - parameter a in equation (19) 

It is seen that both empirical and hydrogenic inputs are required for each edge. When the 
transition is to be computed empirically, s, /3, and a have their appropriate values 
and the corresponding hydrogenic inputs are nonzero arbitrary quantities, except PQN, 
which is zero. If the transition is to be calculated according to the hydrogenic approx- 
imation, cards FB4 have their proper values, while SIG is zero and ESS, AL are arbitrary 
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but finite. Subroutine FBCON is then entered, and a nonequilibrium temperature is 
read if desired, that is, if NEQ = 1. The partition function is then computed. This 
partition function is also used in the LINES subroutine; so, if line radiation for a species 
is to be included, subroutine FBCON must be entered to obtain this function. The free- 
bound emission coefficient is then computed by use of equations (18) and (19). Ten 
equally spaced wavelength points are used to define the spectral variation between succes- 
sive edges, so that at the completion of the calculation there are 10 *NEDGE wavelength 
points in WG(L). After the subroutine output is printed, RACOM is called and WG, G 
are added to the existing WF, F and the main program is re-entered to obtain the next 
input. 

Parameters related to the calculation of atomic and ionic line radiance are read 
as follows: 

Card LI: [TYRAD Format (7A10)] 

Card L2: [NEQ, LIMAX, MION, EION, WMOL Format (313, 2E14.6)] 

LIMAX - number of lines included 

MION - particle degree of ionization 

EION - ionization energy of particle, equivalent K 

WMOL - particle molecular weight 

Cards L3: [WLI(I), ELI(I), FLI(I), IG(I), IU(I), I = 1, LIMAX 

Format (3E14.6, 213)] 

WLI(I) - central wavelength of line, pm 
ELI(I) - lower-level energy of transition, cm - * 

FLI(I) - line f-number 

IG(I) - lower-level degeneracy 

IU(I) - upper-level degeneracy 

The line cards L3 must be read in order of increasing wavelength so that overlapping 
lines may be properly accounted for. Subroutine LINES is then entered, and the standard 
option of reading a nonequilibrium temperature is provided (Format E14.6). The 
integrated line radiance is computed as prescribed by equation (25), and its half -width 
is determined from the combined Stark, Doppler, and resonance widths of equations (30), 
(31), and (32). Each line is defined by 13 wavelength points - six on each side of the 
central wavelength. The line radiance is equated to zero at 18 half-widths on either side 
of the line center. Intermediate wavelengths are spaced in increments which vary 
cubically in such a way that, where the line spectral radiance varies rapidly, that is, 
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near the line center, the increments are smaller than those for the line wings. The pro- 
file is Lorentzian (eq. (26)). The total number of points in a particular line calculation 
is 13*LIMAX. The completion of the line radiance calculation normally ends the radiance 
contribution for an atom or ion and, after RACOM has been called, radiative inputs for a 
different species are read unless the species counter indicates that the final calculation 
has been made. 


The species deck for diatomic molecules begins with card 1A, described previously. 
The number of processes, NRAD, does not include the total number of individual electronic 
band systems for that particular species but only indicates whether band-system calcu- 
lations are to be made at all. For example, if the species calculation for NO is to include 
4 electronic band systems, free-free, and infrared vibration- rotation bands, NRAD 
should be 3. The next molecular input cards are: 

Card 2M: [NEQ, NSTAT, NBD Format (313)] 

NSTAT - number of electronic states used in partition function 


NBD - number of electronic band systems 

Cards 3M: [GM(I), T(I), I = 1, NSTAT Format (F4.0, E14.6)] 

GM(I) - electronic state degeneracy 

T(I) - electronic state energy, equivalent K 

Cards 4M: [v(l), VX(I), B(I), ALPH(I), I = 1, NSTAT 

Format (4E14.6)] 

V(I) - vibrational energy for state I, equivalent K 
VX(I) - vibrational constant for state l, equivalent K 
B(I) - rotational constant for state I, equivalent K 


ALPH(I) - vibration-rotation interaction constant for state I, equivalent K 

Card 5M: [TYRAD Format (7A10)] 

Card 6M: [KUM, KLM, IDU, IDL, JRED, FNQ Format (513, E14.6)] 

KUM - number of Franck-Condon factors used for upper state 

KLM - number of Franck-Condon factors used for lower state 

IDU - identifies upper state: IDU equals the I value of the upper state as 
read on cards 3M 

IDL - similar to IDU but identifies lower state of the transition 

JRED - has value 1 for red-shaded band systems, -1 for violet -shaded systems 

FNQ - band electronic oscillator strength 
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Card 7M: 


|TU, VU, VXU, VYU, BU Format (5E14.6)] 

TU - upper- state energy term 

VU - upper-state vibrational energy term 

VXU - upper-state vibrational first anharmonic energy term 

VYU - upper-state vibrational second anharmonic energy term 

BU - upper- state rotational energy term 

These variables often appear in the literature as T g , o> e , u> e x e , oj g y e , and B g , 
respectively, and are often given in units of cm“l. The input data should be in equiva- 
lent K, which are hc/k (=1.4388) times the given quantities in cm - *. 

Card 8M: [TL, VL, VXL, VYL, BL Format (5E14.6)] 

Variables of card 8M are the corresponding lower-state quantities of card 7M, in the 
same units. 

Cards 9M: [FC(KU, KL), KU = 1, KUM; KL = 1, KLM 

Format (7E10.4)]] 

FC(KU, KL) - band-system Franck-Condon factor array 

Card 10M: [WMIN, WMAX Format (2E14.6)] 

WMIN - minimum of wavelength range for which band contributes, pm 

WMAX - maximum of wavelength range of band 

At this time subroutine SMEAR is called and, as before, if NEQ = 1, a nonequilibrium 
case is indicated. For these band systems, however, three nonequilibrium temperatures 
are provided for: 

TTRON - electronic temperature, K 
TROT - rotational temperature, K 
TVIBR - vibrational temperature, K 

When required, these are read on one card (Format 3E14.6). The diatomic partition 
function is computed according to equations (40), (41), and (42). 

In the smeared-line model, only two wavelength points per band are needed to 
define the spectral profile of the system. One of these wavelengths is at the band head, 
and another is located at some other wavelength toward the red or violet, contingent on 
the band shading. In this program, the second wavelength is located at a small increment 
(0.0001 pm) from the adjacent band head. For each band calculation, the number of 
wavelength points is twice the number of Franck-Condon factors. At each of these wave- 
lengths the band radiance is then computed from equations (37), (38), (39), and (3), and the 
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subsequent G and WG arrays are added to the accumulated radiance function in RACOM, 
as before. 

The free-free contribution for the molecule is computed next, and input parameters 
are identical to those described for the atoms; that is, cards FF1 and FF2. 

For molecules having infrared vibration-rotation bands, the input is 

Card IR1: [TYRAD Format (7A10)] 

Card IR2: [RLAM, VO, VI, BE, ALE Format (5E14.6)] 

RLAM - rotational spectroscopic constant, equivalent K 
VO - vibrational energy term, cm - *- 
VI - vibrational anharmonic constant, cm - * 

BE - rotational constant, cm - *- 

ALE - vibration- rotation interaction constant, cm - * 

Card IR3: (NEQ, NV, NW, WMIN, DW Format (314, 2E14.6)] 

NV - number of vibrational transitions included 
NW - number of wavelength points used 

WMIN - maximum wave-number value, cm - -*- (this corresponds to the 
shortest wavelength) 

DW - wave-number increment, cm"* 

Cards IR4: (FNR(I), I = 1, NV Format (7E10.4)] 

FNR(I) - vibrational transition oscillator strengths 

Subroutine INFRA is entered and, if a nonequilibrium calculation is specified, electronic, 
rotational, and vibrational temperatures are required as inputs in the same manner as 
with the electronic band systems. The band frequencies for the fundamental and first 
two harmonics are computed from equations (45) and (46), and equation (47) is used to 
compute the partition function. Then, for each of the NW wavelength points, the band 
spectral radiance is computed from equation (43) and is transferred to RACOM. 

After all contributions from all species have been computed and added to the 
cumulative spectral radiance in RACOM, the output subroutine, PUTTY, is called. Here 
one card defining the overall spectral range is read: 

Card 01: [RMIN, RMAX Format (2E14.6)] 

Because of the myriad combinations of units presently in use in radiative transfer 
calculations, the output section of the program has been placed in this separate subroutine 
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in order that changes may be made in the form of the output with a minimum of confusion. 
The final output consists of tabulated intensity versus wavelength and optional plots for 
the slab configuration, in addition to the intermediate tabular output. 

Appendixes A and B contain a listing of the program and the printed input data for 
all of the species subdecks used in the sample calculations. 

SAMPLE CALCULATIONS 

Several cases are considered for high-temperature gas mixtures, and comparisons 
are made with similar calculations obtained from other previously developed computer 
codes. Thermally heated COg and air have been chosen for these comparisons. They 
cannot be considered as duplications, since some radiating processes are included in one 
or another of the programs being compared and are omitted elsewhere. Also, basic 
radiative property data (for example, f-numbers, Franck-Condon factor arrays, etc.) are 
not consistently the same for all of the programs used. The primary purpose of the 
sample calculations is to show that the program of this paper computes its solutions in 
the desired manner. 

The equilibrium conditions behind a shock wave in pure COg are used in the first 
example. The thermodynamic state of the resultant mixture, along with the contributing 
radiative processes considered in this program (code name SPRAD), is given in table V. 
References 1 and 2 were used to generate spectral radiative profiles for this mixture and 
are compared in figure 6 with the output of the subject program. In this figure, the 
atomic line contributions of references 1 and 2 are not shown in order to avoid confusion. 
It should also be mentioned that more spectral points could be used in the continuum 
spectrum of the program of reference 1, which would result in greater spectral detail 
than is shown in the figure. Even though the result and spectral integrals (which include 
line contributions for all calculations) are in good agreement (within 10 percent), it is 
in some measure coincidental, since there are relatively severe differences in certain 
spectral regions. 

High-temperature equilibrium air at one-amagat density is also used as a represen- 
tative mixture, and the spectral profiles for a 1-cm slab as computed by program SPRAD 
are shown in figure 7. Reference 4 provides similar spectral profiles with the exception 
of line contributions which are not included. However, certain spectral features may be 
compared; in particular, at the lower temperatures (4000 to 6000 K) the NO infrared 
contribution should be essentially a duplicate calculation since the basic data of Breene 
(ref. 33) are used in the present calculation. Reference 4 verifies this conclusion. 

Reference 3 has also been used to compare the spectral integrals of these air 
calculations. Line contributions have been included in reference 3, but effects of self- 
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absorption in-line overlapping were not accounted for. Spectral profiles are not available 
from this source. The spectral integrals of this program are compared with those of 
references 3 and 4 in figure 8. It should be noted that reference 4 has no inclusion of 
line contributions, and calculations do not extend above 9000 K. The SPRAD results on 
figure 8 are the integrated values of the profiles of figure 7. Differences as large as 
factors of 3.5 exist among the various calculations. 

CONCLUDING REMARKS 

The program described herein provides a means of generating spectral radiative 
properties of gas mixtures in a versatile and convenient manner. It is versatile in the 
sense that the various contributing processes are computed in distinct subroutines which 
lend themselves readily to alterations as required by the advancing state of the art. The 
convenience of using this program is primarily derived from the input arrangement which 
allows the user to omit contributions known in advance to be insignificant, thereby elim- 
inating unnecessary computation. The inherent ability to maintain an up-to-date and 
ready-to-use catalog of punched-card basic data is also advantageous. 

Sample calculations illustrate that the program in its present form is functional 
in the desired manner. While comparison calculations do not exactly duplicate calculations 
of others because of differences in radiative properties included, results are in general 
agreement with outputs of other programs. 

It is expected that this program will undergo continuous change and expansion, 
which is completely in accordance with the manner in which it is arranged. Such changes 
should not alter the basic structure of the program, however, but are anticipated to be 
in the nature of added subroutines and improvements in computational models within 
existing subroutines. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
February 5, 1975 
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APPENDIX A 


PROGRAM LISTING 


PR DORA '-1 S PR AO ( INPUT,JUTPUT , TAP E5= INPUT, TAPE 6 = 0U TP UT ,TAPE21) 

COM IN TT, CUN, EL EC 

COMMON/ KADPT/ G I 650) ,*0(650) ,F( 8000) ,WF( 3000 ) ,NFMAX, NGMAX 
CUMMON/I-R BD/GA! 1 5 J , tA( 15) , GL! 15) , PUN! L5) ,EL( 15) ,EU( L5I , 

1 1 S S t l 5 ) , AL ( 1 3 ) * S I G ( 1 5 ) 

COMMON/ MOL EC/GM( 10 ) , T ( 1 0 ) , 3 ( L 0 ) , V ( 1 0 ) , VX ( L 0 ) , ALPH( 10) , TU, VU , VXU , VY 
lU.BU, TL ,VL , VXL ,VYL,BL, FC I 15, 15 ),WMIN, WMAX 
COMMON/ L I NCR/'* L I ( 50 ) , F L I ( 50 ) , F L I ( 50 ) , I G ( 50 ) , E I ON, M I ON, WMOL , IU( 50) 
COMMON/ INRLO/RL AM f VO ,V 1 , 3E, ALE, FNR(35 ) 

DIMENSION TYP A0( 7 ) 

N S P = 0 

VARIABLES IN COMMON BLOCK PTVAR RELATE TO PLOTTING SUBROUTINE 
INFUPL r 
CALL PSEUDO 
R l A 0 1 , NSPEC » T T , LL EC 
l F )R M A T ( I 3 ,2E14.o ) 

PRINT12.TT 

12 FORMAT! 4X.HTT = F8.C) 

? IF ( NS°.L-O.NSPFC)GO T.J 13 
READ 3, SPEC, ISPFC,NRAD,CON 
3 FORMAT ( A6, 2 13, 614.6 ) 

PrINF 11, SPEC, CON 

ll FORMAT! IH1,4X, A6 , 4X , E 1 4 . 6 , 2 AbHPE R CC/////) 

N S P = >1 S P + 1 

IF! I SPEC . FtJ. 1 ) GO TO 30 
ATONIC SPEC IES INPUT 
N=C 

A READ 5 , TYPAD 
PR I \ I f , TYRAQ 

5 f ■ iRM AT ( 7 4 L 0 ) 

N = N E I 

READ 6 , Nl; 0 , in * , W M I N , DE L W , Z Z 

6 FuRMAT! 12, 15,31. 14.6) 

CALL F FCON ( NC Q , Nw , WM I N , DFLW , ZZ , NS P ) 

IF (N.tr J. NR AD) GO TO 2 
ATOMIC FREE BOUND INPUT 
READ 5 , T Y R A 0 
PRINT ! , T YE AD 

7 F0RMAT(////7A10//) 

N = N E 1 

READ h , NED, NS TAT , NFDGE , X MI M, DE LX 
B FORMAT ( 31 3»2E14.6) 

RcAj ? * ( G A! I), FA! I),t = l,NSTAT) 

9 FORMAT! F4.0»E14.6) 

Rt AD 10, ( GL ( I ) ,PUN( I) ,£L( I ) ,FU( I) ,1 = 1 ,NFDGF) 

LO FOR MAT{F4.Q,3L 14.6) 

READ 11 3, ( ESS! I ) , AH I ) , SIG! I ) , 1 = 1 , NEDGF ) 

113 FORMAT! jF 7. 3) 

CALL f iJCJN(i\ltQ, NS TAT, NE DGE , XM IM, DELX , OA ) 

Ih( N.t.v.NRAO)GO TO 2 
C ATOMIC LINES INPUT 

READ 5, TYRAD 
PRINT 7, TYRAD 


28 



APPENDIX A - Continued 


N = N+ l 

REAL) 3 5 , N E Q * L I MA X , M 1 0 N » E 1 0 N , WMOL 
35 FORMAT ( 31 3, 2E14.6) 

REAJ 34, ( WL I ( 1 ) , EL I( I ) , ELI ( I ) , IGll) , IU< I ), 1 = 1, LI MAX) 
34 FURMAT (3614.6, 213) 

CALL L I NE S ( NE'j * L I M AX , Q A ) 

IE CN.EQ.NKAD) GO TO 2 
3 0 M = 0 

C MOLECULAR SPECIES INPUT 

N B= 0 

READ 14,NFU,NSTAT,NB0 

14 F UK MAI (313) 

KEAU 1 3 , { Of’ ( I ) , T ( I ) , I = l , NS TAT ) 

15 FORMAT ( F4.Q»E 14.6) 

KL AO 16 , ( \M I ) , VX(I) ,B( I ) ,ALPHI I ) » 1 = 1 ,NSTAT) 

l o F 'JK M A F { 4 E 1 4 . 6 ) 

20 REAJ 5 , T /RAD 
PRINT 7 , T Y R A D 

K LA J 1 7 , KJF’.AL M , [UU, I UL, JKED.FNQ 
Nb= NR+ L 

1 7 Kj< MAT( 51 3 , E L 4 . 6 1 

REAJ L 3 , TU, VO, VXU,VVU,BU 
HEAD 13, TL , VL , VXL , VYL , BL 
13 F OR 1AT ( 561*+ .6 ) 

RE AU 1 i , ( ( FC ( K'J , KL ) ,KU=1 ,KUM) , KL = 1, KLM) 

10 F uR 4 I ( 7'1 1 0.4) 

KEAU 2 1 , A M I N, a MA X 

21 F JRMAT ( 2 E 1 4 . 6 ) 

C. ILL S 1EAR (Ni:u,\STAT,KUM,KLM, I DU, IDL, JR ED, END) 

I F ( N.L J.NRAD) GU TO 2 
I F ( N -3 .EG. N130) 22, 20 
2? CONTINUE 

C MULtCULAR ERFE FREE INPUT 

REAJ 5 , T Y R A J 
PUNT 7 , T Y 2 AO 
N = N *■ l 

R t AU o ,NE ), N/J, WM IN.OELW, ZZ 
CALL F F CUN ( NE J » v» , W.3 1 N * J t L W » Z Z ) 

IF ( N.EO. nIKAD ) GU TO 2 
C 3JLECULAR I NFK Ak cl) INPUT 

RtAJ 5 , T Y '< A 0 
PRINT 7 , T Y P A i ) 

N = N + 1 

KEAJ 1 3 ,RL AM , VO , VI , BE , ALE 
READ 4 J , N C 0 , N V , NW , WMI N , DW 
40 FORMAT! 314,2614.6) 

KL A J 1 ) , ( FNR( I ) , I =1 , NV ) 

CALL I N FRA ( NfcO , \ V , NW , rt M I N , Dk ) 

IP{ N.EO. IPAD ) GO TO 2 
C FINAL OUTPUT 

13 CALL PUTTY 
STOP 
END 
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APPENDIX A - Continued 


SU8RJUT IlME FFCON(NEQfNWfWMIN f OELWfZZ»NSP) 

COMMON TT f CON t EL EC 

COMMON/ RADPT/GI650) » WG i 650 ) 9 FI 8000 ), WF < 8000 ), NFMAX, NGMAX 
COMMON /PTVAR/ INC 2) »YM< 5) f XM( 5) ,NXM *NYM f I EC t XMIN»XMAX, YMIN, YMAX, 
il SYMf ITAPE 

DATA IN/20H1250 JENHSSO MS419 / 

DATA YM/50H R AD ( WATTS/MI CRON-CC -SR ) / 

DATA X v l / 5 0 H WAVELENGTH! MICRONS) / 

N XM = 3 
NYM=5 
t FC= L 
XMI N = 0. 

XMAX^O. 

Y M I N = 0 • 

Y M A X = 0 • 

I SYM=14 

ITA1>l = SLTAPE21 
I F ( vJLU.EO. 1)1*2 

1 HEAD 3 » T T R ON 

3 FORMAT ( EL 4. 6) 

GO TJ LL 

2 T TRiJ l \i = T T 

11 PRINT LZrTTRON 

12 F iJR M A T ( 4 X8 H T F M P = t F 8 • 0 t IX INK) 

PRINT 13 

13 FORMAT ( //d/lHNf 5X10H*lAVELENGTH f 3XI3HCRUSS S ECT ION f 2X9HABS COEFF, 

1 6XbHKA») I A-NC Ef iiX9HBLACKB00Yf6X8H INTEGRAL) 

CALCULATION OF WAVELENGTHS 

DO 4 L = 1 > NX 

wG( L)=NMIN+(L-1 )*0ELW 

4 CONTINUE 
L =1 

R AU= J. 

ARt A=0. 

0 E L L A = 0 . 

CALCULATION 'JF RADIATION 
b r 11 K M = K A D 

KAD=( 1.63F-31) < c ELEC*CUN*ZZ*EXP (-1438 / .866 / ( WG ( L ) *T TRON) )/( SQRT 
1 I TTKuN ) *NG ( L )* *2 ) 

1 F ( LL .UT . 0.999 ) GU TO 32 

RAJ=AA J* I J.E- I 0*( WG(U * T TRON— 2 0000. ) **2 / ( WG ( L ) **2. 5 ) 
l L • t-l J*T THUN* *2.5) 

32 CONTINUE 

B3L= ( 1. W06381E+04)/( (WG(L )**5)*(EXRl 1 4387. 8 86/ ( WG ( L J *T TRON ) »—!-)> 

A K S = < A 9 / ft 8 L 

S IGFF=AftS/( ELEC*CONJ 

AR2A=A AF A+. 5*DELLA^( TEKM-t-RAD) 

DELL A = rtG( L+l ) - W G ( L ) 

P K IN r J , L , WG( L ) , SIGFF , ARS ,RAO, B8L , AREA 

5 F !_!•< 3 AT ( NX18»F14.8»i>EL4.6) 

If-(NiP.E0.1il5t 7 

15 F(L)=RAD 
WI (L)=WG(L) 

L =L + I 

IHL.GE.NwJGU TO 8 
GO F.) o 

7 G ( L ) = A AO 
L -L + i 

IFlL.Gc.NWJGO TO 9 
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GJ T i 6 
a nv.o = nw-l 
NFMAX=N WO 

CALL l \i» c JPLT ( I bC , NW.)»WF , i , F , 1 f X*1I N, XMAX t YMI N f YMAX,0. ,NXM, XM» NYM f 
LYMtO) 

on n 10 

9 NWD= nIW- i 
NGM AX = M Wt> 

CALL MFUPLT l L EC f NWO,WG, i tGt i fXMINf XMAX»YMINf YMAXt 0. » NXM, XM,NYM , 
1YM, 0) 

CALL XACOM 
10 RETURN 
EiMU 
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S'JilK.JUr INk F 8 CON (NEUrNSTATtF LDGE,XMIM»DELX»QA) 

COMMON T T » C ON, ELEC 

COMMON/ FR3J/G A ( i 5 ) , t A( i 5 J , GL ( 15) ,PQN( 15) , EL 1 1 5) , EU ( 1 5 ) , 

ItSS(lS) ,AL( 15) ,S IG115) 

COMMON' / R A OP T/ii(650)fRGI650) , FI 3000) f WF (8000) , NFMAX , NGMAX 

COM Mi),, / PTVAR/ 1 N(2 ) ,YM( 5 ) , XM(5 ) , NX M, NYM , 1LC t XM IN , XMA X, YM IN , YMAX , 

L I SY M , L T AP E 

OMFMSIOM 2 Sul 15) ,WEJGE(15 ),CGNC( 15) 

IF! NEy.t 3. I ) 1 , ? 

L HEAD it TTRON 

3 FORMAT (El 4.5) 

GO TO 1 I 

2 T TR0N= T T 
LI PRINT 12, T T RON 
L2 F jkMAT l 4X5H TEM P= , F 8 • 0 » IX l UK ) 

CALCJLATION OF PARTITION FUNCTION 
W A = 0 • 

DU 4 1 = 1,. NIST AT 

w4 = gA#-UA( I )*LXP(-EA{ I) / TTRON ) 

4 C J N T I.NiJE 

CALCULATIOM UP RADIATION 

CR *C<=-6.33E-06 

UN :> 1=L,NPDGL 

I F C EU( I ) .Lt .EL ( I ) ) G 1 J TO 369 

ZSQ ( i ) - — CR A CK* PQN ( I ) **?*(EU ( I )-EL ( I ) ) 

htOGF ( I ) = L .43 3 7886E* 04/ ( EU( U -fcL ( N ) 

CuNC ( I ) = ( CON/G, A ) *CL U ) *EXP (-EL < I ) / TTRON) 

5 CONTINUE 
PRINT 

6 LURMAT( // JXIHN, 5 X 10HWA VfcLENGTH , 3X1 3HCRi:SS S EC T ION , 2X9H ABS COtFF, 
lbXBHRAN I ANCE ,3 X9HBL ACKBUDY , 6 X3H INTEGRAL ) 

I SI T= l 
ARC 4=0. 

DLL L A = 0 . 

RX = 0 • 

XPO=0. 

Uu 10 L-IfNfcDGE 
DO 9 K*L f 10 
POX=XP 3 

X P 0 = X M I K-L) * ( WEDGE ( L ) -XM I M ) /l • 

A d S X= 0 . 

DO 7 I = L , N E 0 G E 

ABSX=ARSX + CONC (I) *S IG(I)*i . F-18* < AL < I ) * ( XPO/ W EDGE ( I ) )**ESS ( I) + 

1(1. -AL <11 )*< XPO/ WEDGE ( U)**<ESS( I) +1 . ) ) +4. 2E-2 5*XPO**3*EU ( I ) **2/ 
LP-J.U I ) * * 3 * C U N C ( I ) 

7 CUNT INJE 
TERM = R:< 

dX= ( L . 19363 3 IP-16 ) / ( ( X PU**5 * l . E- 2<J ) * ( E XP ( 1 . 4 3d 7 8 86F + 04/ { X PO*T T RON ) 
L J - 1 • ) ) 

R X= At3$X * -3X 
SIGFG=4 JSX/CUN 

ARlA=AREA + . 5- DEL LA-(TERM+RX ) 

0LLLA=XP0-Pnx 

PRINT S, I SIT ,XPi)» SIGFBt ABSXt kX,BX » AREA 
6 FORMAT ( 4XI6 f FI4.d f 5EI4.6) 

WG( I SI I* )=XPU 
G ( I SIT ) =RX 
NwD=ISIT 
ISIT=iSIT+l 
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9 CONTINUE 

XMIM=wEf)GE(L)+OELX 
10 CONTINUE 
NGM AX=NwO 

CALL INFOPLT (I EC,NW0,WG»1 ,G, 1 , XM I N, XM AX , YMI N . YMAX , 0. ,NXM, XM,NYM, 
1YM, 0) 

CALL RAC CM 
359 RETURN 
END 
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SUBRUUr INE SMEAR I NE Q , N ST AT , KUM ,KL M , I DU, I DL, JRED,FNQ) 

COMMON T T , CON, ELEC 

COMMON/ RA DPT/ G i 650) ,WG(650),F( 8000 ), WF ( 8000 > ,NFMAX» NGMAX 
COMMON/PrVAR/ IN(2),YM(5),XM(5) , NX M , N YM , l E C , XM I N , XMA X f YM I N, YMAX, 

1 I SYM, IT APE 

COM MU \J / M3 L E C / 3 M (10), ni0),B(10),VU0),VX (LG) ,ALPH( 10) f TU,VU f VXU*VY 
1U,BU,TL,VL, VXL,VYL,BL, FC ( 15,15 ), WMIN, WMAX 
DIMENSION GAM ( 10) , FRAC T ( 10) ,0( 10) 

IF(NEQ.EQ.l) 1,2 

1 READ 3, TfRONfTROT, TV I BR 

3 FORMAT ( 3E 1 A* 6 ) 

GO TO 4 

2 T T RON = TT 
TRDT=TT 

T VI 8R = T T 

4 PRINT 3 , TTRDN, TROT ,TV I BR 

5 FORMAT ( 4X6H TTR ON— , F 8 • 0 , 4X5H TRG T= , F 8 . 0 , 4X 5HT V I 3= , F 8 . 0 ) 

DtLLA=0. 

R A 0 = 0 * 

A R E A- 0 ♦ 

CALCJL AT ION OF PARTITION FUNCTION 
0 1=0. 

DO 9 N= 1 , NS T AT 

G AM ( N ) = ( L • / V ( N ) ) : M2. 3 * t VX(N)+ALPFl{iN)/B(N)-»-8*#8{N)/V{N)) 

Qi N )=GM(N)*tXP < -T (N J/TTRUN) *{ T ROT/BCN) )*( 1 . + GAM ( N ) * T V I SR* TR OT/ { T VI 
1 3R+TRD T ) ) / { l.-EXP(-V(N)/TVIbR)) 

Q T = 0 T + j ( \! ) 

9 CONTINUE 

DO 7 iM= 1 , NS T AT 
FRACUN )=Q(N)/QT 
7 CONTINUE 

CALCULATION OF WAVELENGTHS 
L = 0 

DO 10 KU= 1 9 KUM 
UM=KU-.5 

ENU=TU+ VU*UKM- VXU*UKM* *2+V YU *UKM**3 
DO II <L"I , K L M 
8 KL = KL - . 5 

BNL=rL+VL*B<L-VXL*BKL**2+VYL*BKL**3 
L = L+1 

a G ( L ) = 1 .4 3B78 86F +04/I ENU-ENL ) 

KP=L+1 

WG ( K P ) = WG ( L ) -J RED*. 000 l 
L = L + l 

LI CONTINUE 
10 CUNTINJE 
L M A X = K P 

SORT wAVELENGT HS 
MLX=KP-1 
00 100 N=1,MLX 
J = N 1 

DO LOO K= J , L 1A X 
IF ( WUI N)-aG(K) ) 100,100,200 
2 00 DOL T=WG (N) 

WG ( N) = W G ( K ) 

WG ( K ) =DULT 
LOO CUNTINJE 

CALCULATE ABS COEFF CASE CONSTANTS 
ZEKOL= TL+.5*VL-.25*VXL+. 123*VYL 
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ZfcROU= T J+.5*VU- . 25*VXU+. 12 5*VYU 
SOO=l .43878d6£ +04/ l ZEROU-ZEROL ) 

QVL =0. 

DO 21 KL= l ,53 
6KL = i<L-.3 

T V I B=VL*BKL-VXL*BKL**2+VYL*BKL**3 
I F ( TV 1 1 . LE . 0 . ) GO TO 21 

CVL=wVL+EXP ( -TVI B/TVIBR) 

21 CONTINUE 

8AP=l.273o539t— l 2^ CON* FRACT ( I DL) * SOO*F NU*BL / ( UVL*TKOT ) 

PRINT o? 

5 2 FORMAT! //6XIHN , 5 X 1 OH WA V ELENGTH , 3X L 3HCR OSS SEC T ION , 2X9HABS COEFF, 
16X8HRA0 I ANCE, 5X9 H6L ACKBOOY , 6 Xd H INTEGRAL) 

CALCULATION OF PHI 
L= 1 
NUB= 1 
50 PHI =0. 

DO 15 K'J = L , ,<UM 
Hk.U=XU- . 5 

E VU= VU*3KU-VX J#BKU**2+ VYU*BKU’S c < ! 3 
BVU=BU-ALPH< IDU) * 8KU 
DU 15 M = I,KLM 
BK.L = KL- . 5 

EVL=VL*8KL-VXL*BKL**2+VYL*BKL**3 
8VL = f5L-ALPH( IOL>*BKL 
X=1 . 4 3d 7d86F+04/ WGI L ) 

XVV=TU-TL+F VU-EVL 
IF! JREJ.EO.l) 18,17 

17 IF l X.LT. (XVV-. 0001 ) ) GO TO 16 
GO TO 20 

18 IF! X . G T .( XV V + . OC 0 1 ) ) GO TO 16 

20 CAT=EVL/TV lBR + rWL*(X-XVV) / ( ( BVU-6VL ) *TPOT ) 

IF( CAT. LT.O. » GO TO 16 
IFICAT.GT.100. ) GO TO 16 
U n G=FC( KU,KL ) / ABS( BVU-BVL) 

P H I = PH I + 0<)u*EX Pt -CAT > 

16 CON T I J E 
15 CUNT IN JF 

CALCULATION OF ABSORPTION COEFFICIENT 

AK S = RA p* ! l.-EXPI-l . 4 38 7886 F +04/ I WGI L ) *TTRON ) ) I *PHl /WGI LI 
B iL= 1 . 1906 38 IE- 16/ II WG l L > ** 5 * l . E- 20 I * I EXP I 1 . 4 3B78 86E +04/ 

II WGILI *TT«UNI l-i. )l 
T EK M=RAO 
RAD=ARS*B8L 
SIGMJ=ARS/CUN 

AREA=AREA+.5*DELLA*l TERM+RAD) 

UtLLA= I WG1L + 1) -WGI L) I 

PRINT 40,L,wGt L) , S I GMO , ARS , R AD , BB L , AR E A 
40 FORMAT! 4X 1 d , FI 4 . 8 , 5E 14 . 6 1 

IF ( WGI L) • LT .WMIN.OR.WG IL) .GT .WMAX) GO TO 94 
WGINUB) = WG( L > 

G(NUd) =R AO 
NU8=NUB+L 

94 L = L + 1 

IFIL.GE.LMAXIGO TO 95 
GO T.J 3 0 

95 NW0=NU3-1 

CALL 1NF0PLT II EC,NWD,WG,1,G, l , XM I N, XM AX , YM I N , YMAX , 0. ,NXM,XM,NYM, 
1 Y M , 0 ) 



NGMAX 
CALL 
R E T UK 
END 
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SUBRDUT INF LINFS l NEQ , L I MAX , QA ) 

COMMON TT, CON, ELEC 

COMMON/ L INER/WL I ( SO ) , E L I { 5 0 ) » FL I 1 50 ) , IG ( 50 ) , £ I ON»M ION, WMOL , I U( 50) 
COMMON/ RADPT/GI650) ,WG<650) ,F(8000J , HFI8000) »NFMAX» NGMAX 
COMMON /PT VAR/ I N( 2),YM(51,XM{5) ,NXM, NYM , I EC , XMI N. XMAX, YM I N, YMAX , 
l ISYM, I TAPE 

OIMtNSION FC I 5 0 ) ,HW IDl 50 ) , FUDI 650 ) 

0)131 1 = 1, LI MAX 
FLI [ I ) = 1 . 4388* ELI ( I) 

131 CONTINJE 

C CALCULATION OF k AD I A T ION 

PRINT 6 

6 FOR MAT! / / 8X 1HN , 5 XI OHWA V EL ENGTH, 3X 1 3HCRGSS SECT ION, 2X9HABS CGEFF, 
16X8 HR AO l ANC.F, 5X9 H BLACK BODY, 6X8HINTEGRAL 1 
I F ( Idg.EU.l) GO TO 85 
T T t< ON = T T 
GU TO a 7 
d 5 READ 85 , T TRON 
86 FORMAT l El't.6) 

CuNTINdt 

00 IQ N= 1 , L I MA X 

EC(i>l) = L . J46F-12*CLIN*IG<N)*FLI IN) *EXP(-(EL I ( N 1+ 1438 8/ WL I ( N > ) / TTRON ) 
1/ ( 0 A * V, L I (N 1**3 1 

0 VI = WL I (N » * ( rt I 0N+ 1 ) / ( E I ON- ELI ( N » ) 

H.nUIM =7 .2 3l- 10*ELFC*DV1**2/SQRTITTROM 
1 + 3. 58 IE -07* WL I ( N ) * SURT ( TTRGN/toMJL) 

L + + . 1 563L-22*CLn*FL I (N ) *WL I (K ) **3*S0RT ( IG ( N ) / IU ( N I ) 

1 F ( H H 1 0 ( N ) . G T . 0 . 0 2 ) H W I D I N ) = G . 0 2 
10 CONTINUE 

C CALCULATION OF WAVELENGTHS 

J = 0 

0 0 20 1 =1 , LI MAX 
00 ?l >1 = 1,6 
J-J + l 

WC( J 1 =*LI ( l »-HWIC( I)*( 7-N) **3/ 12 
FUU( J ) = L . 

IF ( N.50. I >FUDI J>=0 
21 CONTINUE 

J = J + I 

rfO ( J ) = .v L I ( I ) 

FUIJI J ) = 1. 

DO 2 3 J = l ,£ 

J=J + 1 

WG( J ) = «L 1 1 I ) + HW IU( I )*N< t *3/L2. 

F UO ( J ) = 1 . 

IFl l.t J • 6 ) FUDIJ ) = 0. 

20 CONTINUE 
M MA X = J 
MLX=M*AX-1 
DO 3 J J = l , ML X 
I=J + 1 

0 1 3U K = 1,0 MAX 

IF (Wi( J).LT.wG(K) 1G0 TO 30 

DDL 1 = <i 3 ( K ) 

DUMB*!. 

W O ( K ) = ri t j ( J ) 

F U I) ( X I =FiJD( J J 
WG( J ) =.) OL T 
FUDIJ ) = DUMB 
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30 CONTINUE 

DO 40 J=1,MMAX 
RAD=0. 

OG 50 1=1, LI MAX 

QWI TN=ULIl I )-18.*HWlD( I ) 

QWI TP=wLI ( I ) +1 8 . *HW ID ( I ) 

IE IWGI .J).LT . Qrt ITN.OR.WGIJ) «GT.QWITP)GO TO 50 

R AD=RAD+EC ( I J * HW I 0 ( I ) / ( 3. 1416* C ( WG( J )-WLI { I I )**2+HWID(I )**2)) 
l*F'JD( J ) 

50 CONTINUE 
G ( J ) =R AD 
40 CONTINUE 
AREA=0. 

DU 9 3 1 = 2* MM AX 
RAD=G( I ) 

XPU= WG ( I ) 

88L = l 1. 190O38E-16) /( ( XPO** 5*1 . E-20) * ( EXP ( l . 43 3 78 86E 4- 04/ ( XPO* TTRON ) 
L ) - 1 - > ) 

AkS=RAJ/83L 
S I GL = AR S/CON 

AP.EA = ARF.A+RAD* (WGU)-WG( 1-1 ) ) 

PRINT 25, I * W G( I ) , SIGL, ARS»RAD»BBL,AREA 
25 FORMAT! 4X I 8 , F l 4. 8 , 5E 14. 6) 

90 CONTINUE 
N rt 0 = M M A X - l 

CALL INF'JPLT (I EC,NVvD,WG,l,G, l , X M I N , XMAX, YM I N , YM AX , 0 . , NXM, XM , NYM , 

1 Y M , 0 ) 

NGMAX=NWD 
CALL RAC OM 
RETURN 
END 
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SUBROUTINE I NlF R A ( N FQ , N V , NW , UM I N, DW ) 

COMMON TT, CON# ELEC 

COMMJN/RAOPT/Gi 650) ,WG (650) ,F( 8000) ,WF<8000 ) , NF MAX, NGMAX 
COMMON/ PT VAR/ I N ( 2 ) , YM C 5 ) , X M ( 5 ) ,NXM, NYM , I EC , XMIN , XMAX, YMIN , YMAX , 

1 I $YM, I r APE 

CUMNJN/ INRLD/RLAM,VO,V 1, 8E,ALE,FNR< 35 ) 

DIMENSION BC ( 3 5 ) ,BROT( 35) v ELV( 35) 

I F ( N E Q . E 0 • i ) LOO, 200 
LOO READ 300, TTRON, TROT, TVIBR 
300 FORMAT ( 3EL4.6) 

GO TO +00 
200 T TR 0N = T T 
T KO T = T T 
T V I 3R- r T 

4 0 0 PRINT 5 00, TTR ON, TROT t TV I BR 

500 FORMAT ( 4XoH T T R UN ■= , F 8 • 0 , 4X5 H TRGT = , F 8 . 0 , 4X5H T V [ B= , F8 . 0 ) 

H= l . A 33 6 

calculation of fundamental band frequencies 

I ~ I 

DO 3 K V ~ 1 , 1 0 

H.Vll)=V3*UV-0.5 >-Vl*(KV-0.5)**2 

BC( I )=VO*(KV+0.5)-VL*d<V+0.5)**2-ELV< N 

BRJTUJ =BL-ALC*l KV-0.5 ) 

1 = 1 + 1 

3 CONTI-mUF 

00 4 nV = 1,10 

EL V ( I ) =V3*l KV-0.5 ) -Vi'- 1 ( KV-0.5 ) **2 
bC(I)=V0*(KV+1.5 >-Vl*< KV + 1. 5)** 2-EL V( l ) 
u R G T { I ) ■= b B- A L F * ( KV-0.5 ) 

1 = 1 + 1 

4 CONTINUE 

DO 5 K V = 1 , 10 

FLVU )=VU*{ KV-0.5 >-Vi*< KV-0.5) **2 

8C( I )=VO*(KV+2.5 )-Vl*(KV+2.5)**2-ELV( I ) 

BROTH) = 3 E-ALL* (KV-0.5) 
l = l+l 

5 CONTINUE 

CALCULATION OF PARTITION FUNCTION 
Q =0 . 

00 63 1=1,1+ 

Q=0+cXP((-{I-.5) *VO*H-Vl*H*< I - . 5 ) ** 2 ) / T V I 8 R > 

50 CONTINUE 

q-o*tkit/rl am 

Al=d.853E-13*C0N/Q 
CALCULATION OF WAVELENGTHS 
DO 11 I = L , N W 
WG { l)=WMIN-(I-l)*DW 
11 CONTINUE 
PRINT 62 

32 FORMAT { //8X1HN, 5 X 1 OH WA VEL E NOTH , 3X 1 3HC R OS S SEC T I ON , 2 X9H ABS CGEFF, 
L6X8HRADI ANCE ,5 X9HBLACK BODY, 6X8 H INTEGRAL) 

CALCULATION OF A 8 SORPT I ON COEFFICIENT 
DiJ 8 I = l , N W 
A 8=0. 

no 7 N = 1 , N V 
0 EL N = k G ( I ) -SC ( N) 

X l = H * E L V < N) /TVIBR 

X2 = H*D EL N** 2/ ( 4 • * BROT ( N) *TROT ) 

Xi=H*0tLN/( 2.* TVIBR) 
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PA=A8$< OFLN) 

Z=X LfXZ+Xi 

IHZ.GT.'JO.) GO TO ? 

Afl=A3f EXP(-Z J*FNK(N)*PA/(4.*BR0T(N)**2) 

7 CONTINUE 

G ( I ) - A 3* ( 1 . - EXP ( -WG ( I ) *H/TTRUN ) I*A1 

8 CONTINUE 
A R F A = 0 . 

W G( L ) = l . E+04/WG ( I) * 

00 9 l = ? » Nw 

WG ( I J = 1 . £ t- 3 4 / w G ( I ) 

ARS=G( i ) 

BF L = (l 19 06. ) /{ { WG( I 1 **5) *( EXP ( 14388. /< wG( I ) *T IRON I » — !.>> 

G ( 1 ) = A 7 S -* 8 P L 
S IGM0=ARS/C0N 

APFA=At<FA + 0.5*(Wb( I )-WG( I-l ) )*(G( I ) +G ( 1-1 ) ) 

PRINT + 0 1 I » W G ( 1 ) ,S iGMO * ARS ,G ( I ) , 8BL, AREA 

40 FORM AT l 4X18 ,H 4. 8,5FL4.6) 

9 CONTI MJF 
NWO=NW- 1 

CALL I IR1PLT ( I tC.NWOfhCf It G, 1, XM I N, XMAX , YM I N , YMAX , 0. , NXM, XM, NYM, 
l Y M » 0 ) 

NGMAX= 1W0 
CALL R A C 0 o 
RETURN 
ENO 


40 



APPENDIX A - Continued 


SUBROUTINE RACOM 

COMMON/ RADP T/G ( 650) « WG ( 650) ,F( 8000) ,WF<8000) t NFMAX » NGMAX 
DIMENSION FF(8000) ,WW(8000) 
l-l 
J = 1 
K = 1 

3 IF(rfF( J).LT.WG(im,2 
2 IF( *G(I) *LE,WF ( J ) )4,5 

5 I F ( WF( J) .LT.WG(I) )6, 2 
l FFC,<) = F(J) 

WW( K) =dF( J) 

IFC J. ED. NFMAX) GO TO 7 
J=J+1 
K =K +• l 
GO TO 3 

4 FF(K)=G( I )+F(J ) — ( W F ( J ) — WG { I))*(F(J)-F(J-1) ) / ( WF { J ) - WF ( J-l ) ) 
W W ( K ) = WG ( I ) 

IFC I • EQ • NGMAX ) GO TO 12 

k=kh 

1=1+1 

GO TJ 2 

6 FF(K) =F ( J) +G(I)- (WG<I)-WF( J) )*(G(I)-GC I-l) )/ (WG(I)-WGC 1-1 ) ) 
« « C K ) = W F ( J ) 

IF( J.EJ.NFMAX)GG TO 10 
K = K+1 
J = J + l 

GO T'J 5 

7 DO 3 1=1, NGMAX 
K = K + 1 

F F ( K ) = G ( l ) 

WW( K) =WG( 1 ) 

8 CONTINUE 
GU TO 0 

10 DO 11 l I - I , NGMAX 
K- K + 1 

FF (K )=G( I I) 

WW(K)=WG(II) 

11 CONTINUF 
GO TO •) 

12 DU 13 JJ=J, NFMAX 
K = «+■ l 

FF(K)=F(JJ) 
w w ( K ) = H F ( J J ) 

13 CUNUNJF 
GO TO ) 

9 K ‘A A X = K. 

DO 14 J = 1 1 K M AX 
f ( J)=FF(J) 

WF ( J J = W* { J ) 

14 CONTINUE 

N F M A X = X M A X 

RETURN 

END 
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SUBROUTINE PUTTY 
COMMON TT fCONj ELEC 

COMMON/ RADPT/G (650) yWG ( 650 ) y F( 8000) ,WF (8000) y NFMAX * NGMAX 
COMMON/ PT VAR/ I N(2 ) y YM ( 5 ) , X M < 5 ) » NX M, NYM , I EC y XM IN , XMAX, YM IN, YMA X , 

1 I SY M f I TAP E 
PRINT 1 

1 FORMAT! lHl f 4XL3H FINAL OUTPUT////) 

PRINT 2 

2 FORMAT! 8X1HN, 5 XL OH*A VE LENGTHy 5X8HR AD l ANC E y 5X9H8L ACK80D Y y 6X8H I NTEGR 
1 AL / ) 

READ 7 » RM IN » RM AX 
7 FORMAT! 2E14,6> 

R A0= 0 • 

AKEA^O. 

OELL A=0 • 

N P T= 0 
I PR I N=0 

00 4 I = 2 » NF MA X 

IF (WF ( I ) .LT .RM IN.UR.WF ( I) •GT.RMAX) GO TO 4 

IF(NPT. EQ • 0 ) NF N= I 

NPT-NPT+I 

BB r JUY=( i. 1906331F— 16) / ( ( wF ( I)**5*L,E-20)*(EXP( I . 4 38 7886E+Q4/ ! W F ( I ) 
L*TT) 1 — i • i ) 

ABSC3=F( I ) / 860 D Y 
TERM=R AD 
R AD=F ( I ) 

DELLA = NF{ I )-WF{ I-l ) 

ARE A=AR EA + . 5*0 ELL A*( TERMER AD) 

F ( I )=ALOGLO! 8B0DY*! I. -EXP {-A3SC0*5. 00) )+L.E-20) 

IF! IPP IN, tQ.O) GO TO 4 
PRINT J V I ,rtF(I) ,F( I J , B80QY y ARE A 
9 FORMAT ! 4X18 y F 1 4 . 8 y 3EI4*6) 

4 CONTINUE 
N trt 0 = N P T 

CALL INFOPLT ( I EC t NWDy wF { NFN ) t I f F ( NFN) y L , XM I N , XMAX , Y1M I N, YM A X y 
10. , NXMyXMyNYMy YM y 0) 

ARE AT =0 . 

NE WN-NF N + I 
F !NF >1) =10,**F( NFN) 

NOMAJ=NFN+NPT- 1 
DU 5 ! - NEW N f NOMAO 
F < I ) = 1 ).**F ! I ) 

AREAT=AKEAT*.5*< wF( I )-WF( 1-1 ) ) *( F { I ) +F < 1-1 ) ) 

5 continue 

PRINT 6 y Ak E AT 

6 FORMAT (4X6HFLUX =yE14.6) 

DO 8 I = NF N y NOM AD 

F m = A LOGIC ( F ( I) *0.80657*WF< I ) **2 ) 
wF { I ) = I # E+04/ ( WF( I ) *8066.) 

8 CONTINUE 

CALL INFOPLTC I EC f NwD, *F < NFN > y l f F < NFN) y I y XM I N , XMA X f YM I N y YM AX y 
I 0. y NXMy XMyNYMy YMyO) 

RET UF N 
END 
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APPENDIX B 


SAMPLE INPUT DATA: SPECIES RADIATIVE PROPERTIES 


C 0 3 3 • 78872 1 E+l 7 

CARBON ATOM FREE-FREE CONTINUUM 
0 51 4 • 990000E— 02 5.000000E-02 2.473000E-02 

CARBON ATOM FREE-BOUND CONTINUUM (HENRY FORMULA) 
0 811 5 . 0 OOOOOE— 02 1.000000E-04 


9# 

4. 

263822E+01 



5* 

1 • 

46673 1 E+04 



1 • 

3. 

1 14902E+04 



5# 

4 . 

S54023E+04 



9 • 

8. 

686884 E+04 



3 . 

8. 

9 1 8371 E+04 



15. 

9. 

221768E+04 



3. 

9. 

907702E+04 



9. 

2. 

000000E+00 

4.263822E+01 

1 .308226E+05 

5. 

2. 

000000E + 0 0 

1 .466731 E + 04 

1 • 308226E+05 

1 • 

2. 

000000E+00 

3. 1 14902E + 04 

1 .308226E+05 

9. 

3. 

00000 OE + 00 

1 . I 62276E + 05 

1 • 30 7560E + 05 

16. 

4. 

000000E+00 

1 .225838E+05 

1 • 30 7560E+05 

25. 

5. 

000000E+00 

1 .255258E+05 

1 • 30 7560E+05 

36. 

6 . 

000000E+00 

1 .271 239E+05 

1 • 30 7560 E + 05 

49. 

7. 

OOOOOOE+OO 

1 .280875E+05 

1 .307660E+05 

64 . 

8. 

OOOOOOE+OO 

1 .2871 29E+05 

1 .307560E + 05 

81 • 

9. 

OOOOOOE+OO 

1 .291 1 47E+05 

1 .307560E+05 

1 00. 

1 • 

00 000 OE + 0 1 

1 .294484E+05 

1 .307B60E + 05 

2. 

000 

3.317 12. 

190 


1 • 

500 

2*789 10. 

300 


1 • 

500 

3.501 9. 

590 


1 • 

000 

1.000 0. 

OOO 


1 • 

000 

1 .000 0. 

OOO 


1 • 

000 

1 • OOO 0 • 

OOO 


1 . 

000 

1 • 000 0 • 

OOO 


1 . 

000 

1.000 0. 

OOO 


1 • 

000 

1 .000 0. 

OOO 


1 • 

000 

1 .000 0. 

OOO 


1 . 

000 

1 .OOO 0. 

OOO 



CARBON ATOM LINE RADIATION 
0 31 0 1.307560E+05 1.201100E+01 


9 

•454400E-02 

3.000000E+01 

2. 700000E-0 1 

9 

3 

1 

•261300E-01 

3.000000E+0 1 

2. 900000E— 02 

9 

9 

1 

. 2 77500E — 0 1 

3.000000E+01 

6 • 30 0000 E— 02 

9 

1 5 

i 

• 280400E-0 1 

3.000000E+01 

2.000000E-02 

9 

9 

1 

• 329300E-0 1 

3.000000E+01 

3.800000E-02 

9 

9 

1 

.431 900E-0 1 

3 . 37350 0E + 04 

1 • 300000E— 0 1 

5 

1 5 

1 

•463330E-01 

1 .01 9400E+04 

9.300000E-02 

5 

7 

1 

.481 770E-0 1 

1 .01 9400E+04 

1 . 1 OOOOOE— 02 

5 

5 

1 

.561 OOOE-O 1 

3.000000E+01 

9. 1 OOOOOE— 02 

9 

1 5 

1 

• 657200E-0 1 

3.000000E+01 

J • 700000E— 0 1 

9 

9 

1 

•751900E-01 

2. 164Q00E+04 

1 .200000E-01 

1 

3 

1 

•930930E-01 

1 .01 9400E+04 

8.200000E-02 

5 

3 

2 

•478560E-0 1 

2. 1 6480 0E + 04 

9.400000E-02 

1 

3 

4 

. 769700E— 0 1 

6.037400E+04 

5.300000E-03 

9 

9 

4 

•932000E-0 1 

6. 198200E+04 

5.500000E-03 

3 

1 
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5.041 700E-0 1 

6.037400E+04 

2.500000E-03 

9 

1 5 

5.052120E— 01 

6. 1 98200E+04 

1 . 1 OOOOOE— 02 

3 

5 

5 • 380240E-'0 1 

6. 1 9820 0E + 04 

7.000000E— 03 

3 

3 

8.335190E— 01 

6. 1 98200E+04 

1 . lOOOOOE-Ol 

3 

1 

9.087600E— 0 1 

6.037400E+04 

3. 1 OOOOOE— 0 1 

9 

9 

9.640600E-0 1 

6.03740CE+04 

1 . OOOOOOE— 0 1 

9 

3 

1 .01 2400E+00 

6.885800E+04 

2 . 620000E— 0 1 

3 

3 

1 .069500E+00 

6.037400E+04 

5. OOOOOOE— 0 1 

9 

15 

1 • 1 33030E+00 

6. 19820OE+04 

4.200000E— 01 

3 

5 

1 • 1 33040E+00 

6.885800E+04 

6 • 300O0OE— 0 1 

3 

5 

1.1641 OOE+OO 

6.972200E+04 

1 . 320000E— 0 1 

15 

1 5 

1 . 1 66400E+00 

7 . 074400E+04 

9.600000E— 01 

3 

9 

1 . 1 88600E+00 

6. 972200E+04 

1 . 400000E— 0 1 

15 

9 

1 • 1 75500E+00 

6.972200E+04 

7. OOOOOOE— 0 1 

15 

21 

1 .2591 OOE+OO 

7. 1 37500E+04 

2.490000E— 01 

9 

9 

1 .698000E+00 

7.2611 OOE+04 

7. 400000E— 01 

5 

7 

C+ 0 3 1 

• 1 4 1 238E+ 1 8 




C+ FREE-FREE 

CONTINUUM 





0 51 5.010000E-02 5.000000E-02 1 .000000E+00 

C+ DUMMY FREE-BOUND (PARTITION FUNCTION DATA) 

0 2 1 5 • 02 1 000E-02 1 .000000E-04 


6. 6.139136E+01 

1 2 • 6.191796E+04 

1* l.OOOOOOE+OO 1 .OOOOOOE+OO O.OOOOOOE+OO 
1.000 1.000 0.000 
CARBON ION LINE RADIATION 
0 20 1 2.829530E+05 1.201100E+01 


6.872500E-02 

4 . 300000E+0 1 

2.600000E-01 

6 

1 0 

8.5841 00E-02 

4.300000E+01 

4.600000E— 02 

6 

2 

9.040900E-02 

4 • 300000E+0 1 

5 • 200000E— 0 1 

6 

6 

1 .01 0200E-0 1 

4 . 30 33 00E+04 

1 .600000E-01 

12 

4 

1 .066000E-0 1 

7 . 49320 0E+ 04 

8. 3 OOOOOE- 02 

1 0 

6 

1 • 036800E-0 1 

4.300000E+01 

5 . 90 000 OE— 02 

6 

2 

1 . 323900E-0 1 

7 . 4932 00E+04 

2.300000E— 01 

10 

1 0 

1 • 335300E— 0 1 

4.300000E+01 

2. 700000E— 01 

6 

10 

2.51 1000E-01 

1 . 1 06520E+05 

1 . 500000E— 0 1 

6 

10 

2 • 837000E-0 1 

9. 64 94 OOE+04 

1 . 300000E— 0 1 

2 

6 

2 . 992600E-0 1 

1 .45551 0E+05 

1 .690000E— 01 

1 0 

14 

3.589300E-01 

1 .81 7400E+05 

1 * 070000E — 0 1 

20 

12 

3 . 920200E-0 1 

1 .31 7320E+05 

1 • 430000E-0 1 

6 

2 

4.267200E-01 

1 .45551 0E+05 

9.400000E— 01 

1 0 

14 

5. 141800E-01 

1 . 670070E+05 

3.420000E— 01 

12 

12 

5.653900E-0 1 

1 . 670070E+05 

1 . 040000E— 0 1 

12 

4 

6.579700E-0 1 

1 . 1 6538 0E + 05 

9.300000E— 01 

2 

6 

6.785600E-0 1 

1 . 670070E+05 

4 . 240000E— 0 1 

12 

20 

7.234400E-0 1 

1 .31 7320E + 05 

5 • 900000 E—0 1 

6 

10 

1 .890200E+00 

1 . 57234 0E + 05 

1 • 1 80000E+00 

2 

6 


C++ 0 3 3.337000E+15 

CARBON DOUBLE ION FREE-FREE CONTINUUM 
0 51 5.025000E-02 5.000000E-02 2. OOOOOOE+OO 

C++ DUMMY FREE-BOUND (PARTITION FUNCTION DATA) 
0 2 1 5.021000E— 02 1.0O0000E-04 

1. O.OOOOOOE+OO 
9. 7.534564E+04 

1. l.OOOOOOE+OO l.OOOOOOE+OO O.OOOOOOE+OO 

1.000 1.000 0.000 
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CARBON DOUBLE ION LINE RADIATION 
0 15 2 5.556850E+05 1.201100E+0I 

9.770260E-02 O.OOOOOOE+OO 8.100000E-01 1 3 

1.175700E-01 5.241900E+04 2.600000E-01 9 9 

2.296B90E-01 1. 02351 0E+05 4.700000E-01 3 5 

3.170160E-01 3.117120E+05 1.470000E-01 1 3 

3.609300E-01 3.177980E+05 3.080000E-0! 9 15 

3.887100E— 01 3.214350E+05 5.800000E-01 15 21 

4.056060E-01 3.242120E+05 5.000000E-01 5 7 

4.122050E-01 3.224030E+05 4.430000E-01 3 5 

4.325700E-01 3.100050E+05 5.000000E-01 3 5 

4.516500E-01 3.177980E+05 1.690000E-01 9 3 

4.648800E-0I 2.382120E+05 7.600000E-01 3 9 

4.662700E-01 3.082830E+05 2.730000E-01 9 9 

5.249600E-01 3.242120E+05 1.280000E-01 5 3 

5 • 696000E-0 1 2.58931 0E+05 4.070000E-01 3 5 

6.740800E-01 3.082830E+05 3.030000E-01 9 15 


C2 1 3 1 .31 0000E+1 3 

0 3 2 

6. 0.000000E+00 

6. 2 . 77776 3E+0 4 

6. 5.766724E+04 

2*361 5557E+03 1 .6790663E+0 1 2. 3489663E+00 2 . 42 1 48 1 2E-02 

2.5728705E+03 2 . 3653685E+0 1 2 • 52 1 7648E+00 2 . 3 1 3572 1 E~02 

1 .5921059E+03 5 . 6486S40E+0 1 1 . 7 1 53238E+00 3 . 48 1 8684E-02 

C2 FOX-HERZBERG BAND SYSTEM 
77311 8. 180000E— 01 

5.766724E+04 1.592106E+03 5.648684E+01 3.100000E+00 1.715324E+00 
O.OOOOOOE+OO 2.361 556E + 03 1.679066E+01 O.OOOOOOE+OO 2.348966E+00 
2000-02 6000-02 1100-01 1700-01 2300-01 2700-01 3100-01 

1400-01 3600-01 5000-02 6700-01 7200-01 7000-01 6500-01 

5300-01 9700-01 1060+00 8700-01 5700-01 3000-01 1100-01 

1240+00 1410+00 7900-01 2000-01 0000+00 1500-01 4200-01 

2020+00 1020+00 3000-02 2900-01 8900-01 0000+00 0000+00 

2410+00 1500-01 6800-01 0000+00 0000+00 0000+00 0000+00 

2130+00 2900-01 0000+00 0000+00 0000+00 0000+00 0000+00 

2.500000E— 01 3.300000E-01 


C2 SWAN BAND SYSTEM 
6 10 2 1-1 1 .500000E-02 

2 • 7777627E+04 2 . 5728705E+03 2 . 3653685E+0 1 -7 . 29034 1 8E-0 1 2. 52 1 7648E+00 

O.OCOOOOOE+OO 2.361 5557E+03 1 .6790663E+0 1 O . OOOOOOOE+OO 2.3489663E+00 
7352+00 2396+00 2445-01 7100-03 0000+00 0000+00 2130+00 

3636+00 3614+00 5983-01 2180-02 0000+00 4312-01 2784+00 

1620+00 4140+00 9832-01 4300-02 7370-02 9158-01 2706+00 

6010-01 4305+00 1274+00 1130-02 2168-01 1293+00 2290+00 

1099-01 5615+0C 1600-03 4220-02 3872-01 1363+00 8880-01 
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0000+00 2000-04 ' 6300-03 6390-02 1216-01 

0000+00 0000+00 3170-02 0000+00 0000+00 

2710-02 0000+00 0000+00 0000+00 0000+00 

0000+00 0000+00 0000+00 0000+00 

4.000000E-01 6.500000E— 01 

C2 MOLECULE FREE-FREE RADIATION 
0 51 5.030000E-02 5.000000E-02 4.946000E-02 

CH 1 3 3.512151E+14 

0 4 2 

2. O.OOOOOOE+OO 
2. 3.330822E+04 

1. 3.733542E+04 

1. 4.578405E+04 

4 • 1 1 7270E+03 9.251484E+01 2.080073E+01 7.683192E-01 

4 • 202735E+03 1.300675E+02 2.145539E+01 9.639960E-01 

3 • 658 1 49E + 0 3 5.378234E + 02 1.854182E + 01 6.978180E-01 

4.063315E+03 1.522250E+02 2.104821E+01 1.070467E+00 

CH BAND (A-X) 

772 1-1 5.000000E-03 

3 * 330822E + 04 4.202735E + 03 1 . 300675E + 02 O.OOOOOOE+OO 2.145539E+01 

O.OOOOOOE+OO 4.11 7270E+03 9.251484E+01 O.OOOOOOE+OO 2.080073E+01 

9996+00 2000-03 0000+00 0000+00 0000+00 0000+00 0000+00 

2000-03 9986+00 1000-03 0000+00 0000+00 0000+00 0000+00 

0000+00 1000-03 9947+00 3200-02 0000+00 0000+00 1000-03 

0000+00 0000+00 2800-02 9819+00 1460-01 0000+00 2300-02 

0000+00 1000+03 0000+00 1310-01 9253+00 4960-01 1280-01 

0000+00 0000+00 0000+00 1600-02 4210-01 8862+00 0000+00 

OOOO+CO 0000+00 0000+00 1000-03 0000+00 0000+00 0000+00 

3.200000E-01 6.200000E— 01 

CH BAND (B-X) 

77311 3.000000E-03 

3 . 733542E+04 3.658149E+03 5.378234E+02 O.OOOOOOE+OO 1.854182E+01 

0. 000000.+OO 4. 1 1 7270E+03 9.251484E+01 O.OOOOOOE+OO 2.080073E+01 

5937+00 1300-02 1845+00 0000+00 8110-01 7000-03 3830-01 

6000-03 2176+00 1660-01 2174+00 1650-01 1610+00 1020-01 

1691+00 1670-01 4000-03 1200-02 5760-01 2150-01 8130-01 

5250-01 1116+00 1620-01 7450-01 1420-01 1000-02 5130-01 

8000-02 1397+00 1800-02 8700-02 3890-01 0000+00 0000+00 

5570-01 6800-02 4530-01 9200-02 4900-02 0000+00 0000+00 

3310-01 4570-01 2440-01 4280-01 3900-02 0000+00 0000+00 

3.000000E— 01 7.000000E— 01 

CH MOLECULE FREE-FREE RADIATION 
0 25 5. 100000E— 02 1.800000E-01 9.956254E-03 

CN 1 3 l.OOOOOOE+OO 

0 3 2 

1. O.OOOOOOE+OO 

2. 1.329690E+04 

1. 3.705169E+04 

2.976453E+03 1.891159E+01 2.733144E+00 2.496318E-02 

2.61 0602E+03 1.853606E+01 2.469700E+00 2.512145E-02 

3.113750E+03 2.91 3570E+0 1 2.834580E+00 3.186942E-02 

CN RED BAND SYSTEM 
77211 7.000000E— 03 

1.329690E+04 2.610602E+03 1.853606E+01 O.OOOOOOE+OO 2.469700E+00 

O.COOOOOE+OO 2.976453E+03 1.891159E+01 O.OOOOOOE+OO 2.733144E+00 

4990+00 3210+00 1260+00 4000-01 1100-01 3000-02 0000+00 


0000+00 0000+00 
0000+00 1000-04 

3800-03 0000+00 
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3700+00 4600-01 2420+00 1950+00 9400-01 3600-01 0000+00 

1110+00 3480+00 1200-01 1000+00 1840+00 0000+00 0000+00 

1800-01 2230+00 2090+00 8900-01 0000+00 0000+00 0000+00 

2000-02 5500-01 2880+00 0000+00 0000+00 0000+00 0000+00 

0000+00 7000-02 0000+00 0000+00 0000+00 0000+00 0000+00 

0000+00 0000+00 0000+00 0000+00 0000+00 0000+00 0000+00 

5.000000E -o 1 2.5Q0000E+00 

CN VIOLET BAND SYSTEM 
7731-1 3.500000E-02 

3*7051 69E+04 3*113750E+03 2*913570E+0l O.OOOOOOE+OO 2.834580E+00 
0 * 0 000 OOE + O 0 2 • 9764 53E + 03 1.891159E + 01 O.OOOOOOE+OO 2.733144E + 00 
9200+00 7900+00 1000-02 0000+00 0000+00 0000+00 0000+00 

7300-01 7870+00 1370+00 3000-02 0000+00 0000+00 0000+00 

5000-02 1210+00 6910+00 1800+00 3000-02 0000+00 0000+00 

0000+00 1400-01 1470+00 6250+00 2110+00 3000-02 0000+00 

0000+00 0000+00 2400-01 1630+00 5700+00 2080+00 2500-01 

0000+00 0000+00 0000+00 2900-01 1710+00 5290+00 2320+00 

0000+00 0000+00 0000+00 0000+00 4500-01 2040+00 4900+00 

3* 400000E-0 1 4 .900 00 0E-0 1 

CN MOLECULE FREE-FREE RADIATION 

0 51 5 • 055000E— 02 1.000000E-01 4.270000E-02 

CO 1 5 l.OOOOOOE+OO 

0 7 3 

1* O.OOOOOOE+OO 
6. 7.0051 02E+04 

3. 8.042972E+04 

6. 8.963567E+04 

2. 9.36288QE+04 

3. 1.205762E+05 

1. 1.250998E+05 

3. 1224734E+03 1 . 9367533E+0 1 2 • 7787324 E+00 2 • 5 1 50 025E-02 

2*50241 31 E+0 3 2.081 92 71E+01 2 • 4 1 86036E+00 2 • 7768620 E- 02 

1 .7524445E+03 1 • 3668492E+ 0 1 1 . 9 1 50276E+00 2 • 30206 1 8E“02 

1 • 6370 393E + 0 3 1 .0969324 E + 0 1 1 .8 150 318E + 00 2 • 44 5Q4 06E-02 

2*1 80 64 24 E + 0 3 2 • 4 8 1 9823E + 0 1 2 . 3 1 875 1 7E + 00 3 . 2 070598E-02 

3* 1624573E+03 0 . 000 0000E+ 00 2 . 9854863E+ 00 4 . 7480024 E-02 
2 • 9956586E + 0 3 O.OOOOOOOE + OO 2.821 4644 E + 00 3 • 8847292E-02 

CO ANGSTROM BAND SYSTEM 
7775-1 3.000000E-02 

1 .250998E+05 2.995659E+03 O.OOOOOOE+OO O.OOOOOOE+OO 2.821464E+00 

9 . 362888E+04 2.180642E+03 2.481982E+01 O.OOOOOOE+OO 2.318752E+00 

8898-01 2505+00 3093+00 2190+00 9767-01 2842-01 5400-02 

1816+00 1757+00 8330-02 8323-01 2320+00 2020+00 9010-01 

2106+00 3039-01 7102-01 1232+00 1870-02 1052+00 2280+00 

1834+00 4200-02 1173+00 2690-02 9680-01 6931-01 1153-01 

1340+00 5214-01 5215-01 4339-01 6583-01 1591-01 1123+00 

8706-01 9553-01 2650-02 8487-01 2500-03 8492-01 3957-01 

6211-01 1066+00 1132-01 5222-01 3376-01 3583-01 6574-01 

2.200000E-01 1 .OOOOOOE+OO 

CO HOPFIELD-BIRGE (B-X) BAND SYSTEM 
2961-1 1.200000E-02 

1.205762E+05 3.162457E+03 O.OOOOOOE+OO O.OOOOOOE+OO 2.985486E+00 

O.COOOOOE+OQ 3. 122473E+03 1 . 936753E+0 1 4.431504E-02 2.778732E+00 

6468+00 2836+00 2597+00 1632+00 6832-01 2998+00 1767-01 

1552+00 5070-02 6032-01 1580-02 2289-01 5200-03 8930-02 

1800-03 3550-02 6000-04 1430-02 

1 • 1 50000E-0 1 1 .600000E-01 
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CO FOURTH positive band system 
1414 5 1 1 1 .500000E— 01 

9 • 3628880E + 0 4 2 . 1 806424E+03 2 . 48 1 9B23E + 0 1 0. OOOOOOOE+OO 2. 3 1 875 1 7E+00 

O.OOOOOOOE+OO 3.1 22 4 734 E+ 03 1 . 9367533E+0 1 4 . 43 1 4689E-02 2 . 7787324E+00 


1 1 32 + 00 

2161+00 

2300+00 

1813+00 

1 188+00 

6885-01 

3668-01 

1843-01 

8897-02 

4184-02 

1 937-02 

8894-03 

4079-03 

1877-03 

2609+00 

1 549+00 

1218-01 

2049-01 

8728-01 

1 230+00 

1 160+00 

8806-01 

5837-01 

3537-01 

2016-01 

1 1 02-01 

5864-02 

3067-02 

2848+00 

3051-02 

90 1 3-01 

1 1 70+00 

3443-01 

3236-03 

3211-01 

7347-01 

9097-0 1 

8453-01 

6633-01 

4664-01 

3045-01 

1 888-0 1 

1 963+00 

7644-0 1 

1 161+00 

6454-03 

5763-01 

9125-01 

4220-0 1 

2234-02 

9899-02 

4063-01 

6476-01 

71 75-0 1 

6509-01 

5204-01 

9604-01 

1 931+00 

5083-02 

8957-0 1 

6656-01 

1 384-04 

4290-01 

7389-01 

4462-0 1 

7969-02 

1462-02 

1948-01 

4177-01 

5531-01 

3554-01 

1 857+00 

5718-01 

8416-0 1 

6171-02 

8236-01 

4226-01 

4267-04 

3278-01 

6081-01 

4505-01 

1450-01 

1 71 4-03 

6720-02 

1034-01 

1 083 + 00 

1 650+0 C 

4713-03 

9803-01 

2268-01 

2178-0 1 

6988-01 

2974-0 1 

3829-04 

2418-01 

4970-01 

441 8-01 

2100-01 

2428-02 

4465-0 1 

1 668+00 

6766-01 

5141-01 

3155-01 

7288-0 1 

2902-02 

3061-0 1 

5834-01 

2373-01 

1 1 03-04 

1630-01 

3919-0 1 

4686-03 

1380-01 

9963-01 

1588+00 

2720-02 

9438-01 

5930-03 

5818-01 

4267-01 

1 669-04 

3252-01 

4944-0 1 

2159-01 

4229-03 

7526-04 

3374-02 

4131-01 

1 486+00 

891 1-01 

21 09-01 

6352-0 1 

3750-01 

9839-02 

5888-01 

2374-01 

1 094-02 

3020-01 

4288-01 

1015-04 

6664-03 

1283-01 

8454-0 1 

1 568+00 

1 802-01 

7283-01 

1 01 0-01 

6822— C 1 

5967-02 

2480-0 1 

4991-01 

1 438-01 

1 718-02 

1 1 57-05 

1 081-03 

31 1 5-02 

3372-01 

1291 +00 

1 1 30 + 00 

2433-02 

8281-01 

6076-02 

4350-0 1 

4170-01 

5751-05 

3144-01 

4100-01 

1 1 1 8-06 

1 458-04 

6072-03 

101 1-01 

676 1 -01 

151 0+00 

4547-01 

4048-01 

4075-0 1 

4454-01 

4708-02 

5187-01 

2062-01 

2127-02 

9187-08 

1646-05 

9671-04 

2372-02 

2528-01 

1 078+00 

1321+00 

2740-02 

7797-01 

2455-02 

6218-0 1 

8492-02 

2073-01 

4459-0 1 


1.200000E-01 2.400000E-01 

CO MOLECULE FREE-FREE CONTINUUM 


0 51 5.000000E— 02 5.000000E-02 4.000000E-02 

CO INFRARED vibration-rotation BANDS 
2.779000E+00 2.169560E+03 1.345300E+0I 1.930240E+00 1 .747000E-02 

0 30 300 1 .OOOOOOE+C4 2.866666E+01 


1 090-04 

2210-04 

3360-04 

4560-04 

5800-04 

7100-04 

8450-04 

9850-04 

1 1 30-03 

1270-03 

7500-07 

2230-06 

4370-06 

7210-06 

1 070-05 

1500-05 

2000-05 

2560-05 

3180-05 

3850-05 

2820-09 

1250-08 

8710-07 

3440-08 
1 1 80-06 

7550-08 

1 450-07 

2520-07 

4050-07 

61 1 0-07 


C0+ 1 3 l.OOOOOOE+OO 

0 3 2 

2. O.OOOOOOE+OO 
4. 2 • 983 1 60E+04 

2. 6.600920E+04 

3.47371 00E+03 2 . 1 8 1 9000E+0 1 2 . 8449000E+00 4 . 3520000E-02 
2.2475500E+03 1 . 9470000E+0 1 2 .2869000E+00 2 . 7940000E-02 

2.49521 00E+03 4 . 0 1 8200 OE+O 1 2 . 5898000E+00 2 . 7280000E-02 

CO+ first NEGATIVE BAND SYSTEM 
10 10 3 1 1 1.700000E— 02 

6.6009200E+04 2 . 4952 1 00E+03 4 .0 1 82000E+0 1 4 . 7240000E-0 1 2 . 58980 OOE+OO 

O.OOOOOOOE+OO 3.47371 0OE+03 2 . 1 8 1 9000E+0 1 - 1 . 0072000E-03 2 . 8449000E+00 
5315+00 3125+00 1115+00 3247-01 8761-02 2351-01 6548-03 

1938-03 6161-04 2106-04 3394+00 5651-01. 2578+00 1984+00 
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9334-01 

3557-01 

1242-01 

4255-02 

1486-02 

5412-03 

1 039+00 

31 80+00 

1763-01 

9127-01 

1864+00 

1445+00 

7765-01 

3518-01 

1478-01 

6087-02 

2142-01 

2170+00 

1464+00 

1230+00 

2037-02 

9643-01 

1438+00 

1 1 34+00 

6772-01 

3522-01 

3383-02 

7490-01 

2550+00 

1 728-0 1 

1646+00 

3359-0 1 

1413-01 

8642-01 

1 139+00 

9377-01 

4328-03 

1748-01 

1471+00 

1927+00 

1 306-01 

9965-0 1 

1 028+00 

7692-02 

1929-01 

7072-01 

4649-04 

3087-02 

5006-01 

2009+00 

8076-0 1 

8407-01 

1652-01 

1 102 + 00 

6058-01 

2128-02 

4280-05 

4378-03 

1 203-01 

101 1+00 

1 977+00 

5947-02 

1233+00 

7172-02 

4898-0 1 

9238-01 

3420-06 

5162-04 

2228-02 

331 1-01 

1548+00 

1306+00 

1665-01 

8166-01 

6065-01 

1 171-02 

2395-07 

5168-05 

1476-01 

3333-03 

8754-01 

7975-02 

7009-0 1 

181 1 +00 

4288-01 

7613-01 


1 •700000E-01 3.500000E-01 

CO+ COMET TAIL BAND SYSTEM 
10 10 2 1 1 7.240000E-03 

2 • 983 1600E+Q4 2 • 247550 0E+03 1 . 9470000E+0 1 1 . 8S00000E-02 2 • 2869 000E+0 0 

O.OOOOOOOE+OO 3.47371 00E+03 2 • 1 8 1 900QE+0 1 - 1 • 0072000E-03 2 • 84490 00E+0 0 


4274-01 

1 151+00 

1698+00 

1 824+00 

1600+00 

1219+00 

8369-0 1 

5316-01 

3182-01 

1818-01 

1 520+00 

1932+00 

9746-0 1 

1272-0 1 

5796-02 

4537-0 1 

8359-01 

9854-0 1 

9185-01 

7381-01 

2499+00 

8046-01 

3228-02 

7483-01 

9618-01 

4719-01 

4973-02 

5155-02 

3212-01 

5927-0 1 

2516+00 

5677-03 

1 083 + 00 

6884-01 

5787-03 

3336-01 

71 97-01 

5649-01 

1 880-0 1 

2345-03 

1 734+00 

9166-0 1 

8454-0 1 

4186-02 

7785-01 

5594-0 1 

2730-02 

1681-01 

5238-01 

5427-01 

8654-0 1 

1 905+00 

3505-03 

9742-01 

3532-0 1 

9275-02 

6355-01 

4599-01 

4066-02 

9005-02 

3233-01 

1 760+00 

7225-0 1 

6206-01 

1993-01 

7857-0 1 

1 316-01 

1385-01 

5420-0 1 

3865-0 1 

9207-02 

9992-0 1 

1 724+00 

1 054-02 

9837-01 

5569-02 

4380-0 1 

5538-01 

3702-02 

1 694-01 

2018-02 

3897-0 1 

161 3+00 

9145-0 1 

2964-01 

5383-01 

5073-01 

2289-02 

52 1 7-01 

3607-01 

34 1 4-03 

1097-01 

8830-0 1 

1 724+OC 

1638-01 

8163-01 

4949-02 

691 7-0 1 


1343-01 1486-01 

3.000000E-0 1 8.500000E-01 

CO + FREE-FREE RADIATION 

0 51 5#040000E-02 5.000000E-02 1.000000E+00 

H 0 3 1*51 OOOOE+1 7 

HYDROGEN ATOM FREE-FREE CONTINUUM 

0 50 5*250000E-02 l.OOOOOOE-Ol 7.315000E-02 

HYDROGEN ATOM FREE-BOUND CONTINUUM 

066 5.100000E-02 1.000O00E-04 

2* 0.000000E+00 

8. 1.1 83540E+05 

1 4 . 1 .40271 0E + 05 

22. 1.479430E+05 

32. 1.51 4940E+05 

44. 1.534220E+05 

2. l.OOOOOOE+OO O.COOOOOE+OO 1.578060E+05 
8. 2.000000E+00 1.183540E+05 1.578060E+35 

14. 3.000000E+00 1*402710 E+05 1.578060E+05 

22. 4.000000E+00 l*479430E+05 1.578060E+05 

32. 5.000000E+00 1.514940E+05 1.578060E+05 

44. 6* OOOOOOE-fOO 1.534220E+05 1.578060E+05 

1 *000 1 .000 0.000 
1 *000 1 * 000 0 • 000 
1 *000 1 *000 0.000 
1 .000 1 .000 0.000 
1 .000 1 .000 0.000 
1 • OOO 1 .000 o • 000 
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HYDROGEN ATOM LINE RADIATION 


0 20 0 1.578060E+04 1.008000E+00 



9.378030E— 02 

O.OOOOOOE+OO 

7.799000E— 03 

2 

12 

9.497430E-02 

O.OOOOOOE+OO 

1 . 394000E— 02 

2 

50 

9.925370E— 02 

O.OOOOOOE+OO 

2 • 8990 00 E— 02 

2 

32 

1 . 025720E— 0 1 

O.OOOOOOE+OO 

7.91 0000E-02 

2 

1 8 

1 .215670E— 01 

O.OOOOOOE+OO 

4. 1 62000E— 0 1 

2 

8 

3.970070E— 0 1 

8.225900E+04 

1 .270000E-02 

8 

98 

4. 1 01730E-0 1 

8.225900E+04 

2.209000E— 02 

8 

72 

4.340360E— 01 

8.225900E+04 

4.467000E— 02 

8 

50 

4.861 320E— 0 1 

8.225900E+04 

1 . 193000E-01 

8 

32 

6.562800E — 0 1 

8.225900E+04 

6.40700OE— 01 

8 

1 8 

9.545980E-0 1 

9. 749200E+04 

1 .604000E-02 

18128 

1 .00494 OE + 0 0 

9 . 74920 0E + 04 

2.768000E— 02 

18 

98 

1 .09381 OE+00 

9 . 7492 0 OE + 04 

5.584000E— 02 

18 

72 

1 .28181 OE+00 

9. 749200E+04 

1 . 506000E— 0 1 

18 

50 

1 .81 741 OE+00 

1 . 02824 0E + 05 

1 .870000E— 02 

32 162 

1 .8751 00E+00 

9. 749200F+04 

8.421000E— 01 

18 

32 

1 . 944560E + 0 0 

1 .028240E+05 

3.230000E— 02 

32128 

2. 1 65500E+00 

1 . 02824 0E + 05 

6.549000E— 02 

32 

98 

2.625200E+00 

1 .028240E + 05 

1 . 793000E — 0 1 

32 

72 

4.051200E+00 

1 .02824 0E+ 05 

1 . 038000E + 00 

32 

50 


H+ 0 1 5.260000E+16 

HYDROGEN ION FREE-FREE CONTINUUM 
0 48 5.150000E-02 1.000000E-01 l.OOOOOOE+OO 

H2 1 3 2.357000E+12 

0 3 2 

1 . O.OOOOOOE+OO 

1. 1.319220E+05 

2. 1.439420E+05 

6.323800E+03 1.696300E+02 8.748000E+01 4.306000E+00 

1.952300E+03 2.867800E+01 2.880000E+01 1 . 7 1 6900E+00 

3.514600E+03 9.644300E+01 4.509200E+01 2.33Q500E+0C 

H2 LYMAN BAND SYSTEM 
4 14 2 1 1 5.000000E-01 

1.31 922OE+0 5 1.952300E + 03 2.867800E + 01 5.797000E-0I 2.880000E+01 

O.OOOOOOE+OO 6 . 32 380 0E + 03 1.696300E + 02 4.173000E-01 8.748000E+01 

7240-02 3214-01 7476-01 1211+00 3630-01 1004+00 1318+00 

1013+00 9947-01 1500+00 8068-01 7843-02 1816+00 1056+00 

1740-02 3713-01 2353+00 1467-01 5008-01 7850-01 2186+00 

2471-01 1055+00 7121-02 1423+00 1510+00 2582-01 4748-01 

6144-01 2208+00 3014-01 8990-01 1579-01 1485+00 1960+00 

5371-03 1927-02 4679-01 2177+00 1410+00 5652-04 5213-02 

7925-01 2655+00 4436-06 4187-04 6404-02 9883-01 4491-07 

8971-05 2222-04 3255-02 1060-07 9728-07 1667-04 4856-03 

8.00000OE-O2 1 • 90000 OE — 0 1 

H2 WERNER BAND SYSTEM 
4 14 3 1 1 3.900000E-0 1 

1.439420E+05 3.345400E+03 9.64430CE+01 O.OOOOOOE+OO 4.509200E+01 

O.OOOOOOE+OO 6 . 32380 0E + 03 1.696300E + 02 4.173000E-01 8.748000E+01 

1221+00 1868+00 1790+00 1403+00 3364+00 1434+00 1266-01 

6821-02 3592+00 3669-02 1163+00 1066+00 1901+00 2324+00 

8277-01 3844-02 5246-01 3315+00 4172-01 1445+00 7144-02 

1564+00 3223+00 4631-02 3872-03 2941-01 2723+00 2157+00 

3554-05 1793-02 6903-01 3628+00 5523-07 8195-05 4344-02 

1210+00 2138-09 8631-06 8292-06 6826-02 1009-09 2946-08 
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4428-06 4679-04 9523-11 1353-08 2403-06 6928-05 4728-11 

5313-09 1579-08 2207-05 5023-11 1265-09 1136-08 3221-06 

8.000000E-02 2.000000E— 01 

HYDROGEN MOLECULE FREE— FREE CONTINUUM 
0 48 5.075000E-O2 1.000000E-01 1.070000E-02 

HE 0 3 1 .890000E+1 9 

HELIUM ATOM FREE-FREE CONTINUUM 
0 50 5.000000E-02 1.000000E-01 2.256345F-02 

HELIUM ATOM FREE-BOUND CONTINUUM 
032 5.500000E-02 1.000000E-04 

1. O.OOOOOOE+OO 
3. 2 • 29992 OE+O 5 

1. 2.392320E+05 

3. 2.000000E+00 2.299920E+05 2.853210E+05 

1. 2.000000E+00 2.392320E+05 2.85321 0E+05 

1 .000 1 .000 0.000 
1 .000 1 .000 0.000 
HFL IUM ATOM LINE RADIATION 
0 41 0 2.85321 0E+04 4.003000E+00 

5.370300E-02 O.OOOOOOE+OO 7.340000E-02 1 3 

5.843340E-02 O.OOOOOOE+OO 2.762000E-01 1 3 

3.447590E— 01 1.662780E+05 1.280000E-02 1 3 

3.61 3640E-0 1 1.662780E+05 2.210000E-02 1 3 

3 . 964730E-0 1 1.662780E+05 5.070000E-02 1 3 

4.437550E-01 1.711350E+05 3.080000E-03 3 1 

4 .47 1 500E-0 1 1.690870E+05 1.250000E-01 9 15 

4.921930E— 01 1.711350E+05 1.220000E-01 3 5 

5.01 5680E-0 1 1.662780E+05 1.514000E-01 1 3 

5.047740E— 0 1 1.711 350E+05 8.340000E-03 3 1 

5.87570OE— 01 1.690870E+05 6.090030E-01 9 15 

6.6781 50E-0 1 1.711350E+05 7.110000E-01 3 5 

7.065300E— 01 1.690870E+05 6.930000E-02 9 3 

7.281350E— 01 1.711350E+05 4.800000E-02 3 1 

1.083000E+00 1.598560E+05 5.391000E-01 3 9 

1 .10131 OE+OO 1.848650E+05 5.210000E-02 1 3 

1 . 1 9691 OE+OO 1.855650E+05 1.230000E-01 9 15 

1.278500E+00 1.861020E+05 1.580000E-01 15 21 

1 .279030E+00 1.861050E+05 1.580000E-01 5 7 

1.296840E+00 1 . 862 1 00E+05 1.390000E-01 3 5 

1.508370E+00 1.848650E+05 1.400000E-01 1 3 

1.700200E+00 1.855650E+05 4.820000E-01 9 15 

1.868600E+00 1.861020E+05 1.020000E+00 15 21 

1.869690E+00 1.861050E+05 1.010000E+00 5 7 

1.908940E+00 1 • 862 1 00E+05 6.470000E-01 3 5 

2.058130E+00 1.662780E+05 3.764000E-01 1 3 

2.11 2000E+00 1.855650E+05 1.450000E-01 9 3 

2.11 320 OE+OO 1 . 862 1 00E+05 1.030000E-01 3 1 

2.161 700E+00 1.914470E+05 7.190000E-02 5 7 

2.472700E+00 1.912170E+05 1.210000E-01 9 15 

2.618500E+00 1.914450E+05 1.870000E-01 15 21 

2.619800E+00 1.914470E+05 1.870000E-01 5 7 

2.6531 OOE+OO 1.914930E+05 1.520000E-01 3 5 

3.329900E+00 1.909400E+05 1.510000E-01 1 3 

3.702600E+00 1.912170E+05 4.420000E-01 9 15 

4 . 036500E+00 1.914450E+05 8.880000E-01 15 21 

4 . 039600E+00 1.914470E+05 8.870000E-01 5 7 

4.121600E+00 1.914930E+05 6.490000E-01 3 5 
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4 • 294 VOOE+OO 1.832370E+05 8.960000E-01 3 9 

4.605300E+00 1.914930E+05 1.590000E-01 3 1 

4 • 693600E+00 1.912170E+05 2.230000E-01 9 3 

HE+ 0 1 1.390000E+18 

HELIUM ION FREE-FREE RADIATION 
0 48 5.150000E-02 l.OOOOOOE-01 1.000000E+00 

HG 0 3 6.440700E+12 

HG F-F 

0 51 1.990000E-01 1.000000E-02 2.000000E-02 

HG F-B 

055 2.000000E-01 1.000Q00E-04 

1. O.OOOOOOE+OO 

9. 6.012500E+04 

2. 7.779400E+04 

2. 8.970900E+04 

1. 9.198000E+04 

9. 6.000000E+00 6.012500E+04 1 .2112408+05 

3. 6.000000E+00 7.779400E+04 1.211240E+05 

3. 7.0C0000E+00 8.97090OE+04 1.211240E+05 

1. 7.000000E+00 9 • 1 980 00 E +04 1.211240E+05 

128. 8.000000E+00 1.192310E+05 1.211240E+05 


1 .000 

1 .000 

0.000 

’1 .000 

1 .ooo 

0.000 

1 .000 

1 .000 

0.000 

1 .000 

1 .ooo 

0.000 

1 .000 

1 .000 

0.000 


HG LINFS 

0 19 0 1.211 240E+05 2.006100F+02 

2.536520E-01 O.OOOOOOE+OO 3.400000E-01 1 3 

2 • 752780E— 0 1 3.764500E+04 1.200000E-01 1 3 

2 . 893600E— 0 1 3.941200E+04 7.666667E-02 3 3 

2.967280E-01 3.764500E+04 2.900000E+00 1 3 

3 . 02 1 500E-0 1 4.404300E+04 1.080000E-01 5 7 

3. 1 25660E— 0 1 3.941200E+04 3.666667E-02 3 5 

3.131 550E-0 1 3. 94 1200E+04 3.066667E-01 3 3 

3.131830E-01 3.941200E+04 3.08O000E-01 3 5 

3.341480E-01 4 . 404300E+04 4.040000E-02 5 3 

3 . 650 1 50E— 0 1 4 . 404300E+04 2.600000E+00 5 7 

3.654830E— 01 4.404300E+04 7.000000E-01 5 5 

3.662880E-01 4.404300E+04 7.400000E-02 5 3 

3.663280E-01 4.404300E+04 2.020000E-01 5 5 

4 . 046560E-0 1 3.764500E+04 8.800000E+00 1 3 

4.077810E— 0 1 3,94 1200E+04 2.133333E-01 3 1 

4.358350E— 0 1 3.941200E+04 8.000000E-01 3 3 

5.460740E-01 4.404300E+04 7.600000E+00 5 3 

5.769590E-01 5.406900E+04 1.433333E+00 3 5 

5.790650E-01 5.406900E+04 1 . 700000E+00 3 5 

N 0 3 1 .892000E+ 1 9 

NITROGEN ATOM FREE -FREE RADIATION 
0 51 4.990000E— 02 1.000000E-01 1.796849E-02 

3 NITROGEN ATOM FREE-BOUND RADIATION 
0 3 14 5. OOOOOOE— 02 1 .000000E-04 

4. O.OOOOOOE+OO 
10. 2.766265E+04 

6. 4. 149500E+04 
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6 

4 

10 

6 

10 

6 

9 

16 

35 

36 
49 
64 
81 

1 00 


а. ooooooE+oo 
2 • OOOOOOE+OO 
2.000000E+00 
2.000000E+00 
2* OOOOOOE + OO 
2 • OOOOOOE+OO 
3 • OOOOOOE+OO 
4. OOOOOOE+OO 
5# OOOOOOE+OO 

б. 000000E+00 
7. OOOOOOE+OO 
8 • OOOOOOE+OO 
9 • OOOOOOE+OO 
1 •000000E+01 


4*1 49500E+04 
0 •OOOOOOE+OO 
2.766200E+04 
4# 1 49500E+04 
2 • 7662 00 E+ 04 
4 « 1 49500E+04 
l • 499089E+05 
1 #581 070E+05 
1 •619016E+05 
1 • 63962 9E+ 05 
1 • 652057E+05 
1 • 660 1 24E+05 
1 .665654E+05 
1 * 6696 1 0E+05 


ooo 

4.727 

2.030 

000 

4.287 

1 1 .420 

000 

4.826 

5.020 

ooo 

5.112 

2.870 

500 

3*847 

4.410 

500 

4.337 

4.200 

OOO 

1 *000 

0.000 

ooo 

1 .000 

0.000 

ooo 

1 .000 

0.000 

ooo 

1 .000 

0.000 

ooo 

1 .000 

0.000 

ooo 

1 .000 

0.000 

ooo 

1 .000 

0.000 

ooo 

1 .000 

0.000 


2. 156775E+05 
1 •687760E+05 
1 » 90683 5E+ 05 
1 .906835E+05 
1 .687760E+05 
1 •687760E+05 
l « 68647SE+05 
1 • 686475E+05 
1 .686475E+05 
1 • 6864 75 E+ 05 
1 •686475E+05 
1 • 6864 7 5 E+ 05 
1 • 686475E+05 
1 • 686475E+05 


NITROGEN ATOM LINE RADIATION 
0 26 O 1 .686470E+05 1 •400800E+01 

1.134600E-01 O. OOOOOOE+OO 1.300000E-01 4 12 

1*1 67900E—0 1 1.922800E+04 3.400000E-02 10 14 

1 . 1 769O0E-0 1 I.922800E+04 1 . 400000E-02 10 6 

1.199900E-01 0*OO00Q0E+00 3.500000E-01 4 12 

1 *2433008-0 1 1.922800E+04 l.lOOOOOE-Ol 10 10 

1.310700E-01 2.884000E+04 5.600000E-02 6 10 

1*31 9500E — 0 1 2 • 884 OOOE + 04 3.400000E-02 6 6 

1.41 1940E-01 2 • 884 00 OE + 04 2.600000E-02 6 10 

1.493300E-01 1.922800E+04 1 .1000008-01 10 6 

1 .743600E-0 1 2 • 884 OOOE + 04 9.100000E-02 6 6 

8.211800E-01 8.333700E+04 2.310000E-02 12 12 

8*61 7500E-0 1 8*61 930 OE + 04 3. 1Q0000E-01 6 6 

8.691600E-01 8 * 33370 0E + 04 3.580000E-02 12 20 

9 • 04760QE— O 1 9.966300E+04 4*670000E-01 10 14 

9.0504008-01 9 . 35820 0E+04 9.450000E-01 2 6 

9.395300E-01 8.619300E+04 4.780000E-01 6 lO 

9.829200E-01 9.483900E+04 1.440000E-01 20 20 

1*01 1 700E + 0 O 9 * 48 390 0E + 04 8.020000E-01 20 28 

1.052S00E+00 9.551 100E+04 6.880000E-01 12 20 

1.059500E+00 1 . 1071 30E+05 7.310000E-01 14 18 


1.070800E+00 9.551 100E+04 2.160000E-01 12 12 
1.129000E+00 9 • 48 390 0E+04 1.680000E-C1 20 12 
1.204300E+00 9 . 68340 0E+04 1.300000E-01 10 10 
1.221700E+00 9 . 55 1 1 0 0E+04 1.71O00OE-01 12 12 
I.235000E + 00 9 . 67520 0E + 04 8.500000E-01 4 12 
1.246900E+00 9.683400E+04 7.100000E-01 10 14 
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N+ 0 3 2.806000E+15 

NITROGEN ION FREE-FREE RADIATION 
0 51 5.020000E-02 1.000000E-01 1 .000000E+00 

NITROGEN ION DUMMY F-B 
0 3 1 5.021000E-02 1 .000000E-04 

9. 1.285000E+02 

5. 2.203600E+04 

1* 4.703000E+04 

0. O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 
0.000 0.000 0.000 
NITROGEN ION LINE RADIATION 
0 39 1 3.435151E+05 1.400800E+01 


9. 1 63400E-02 

8. 9300C0E+0 1 

2 . 2 000 OOE-O 1 

9 

9 

1 .0851 00E-01 

8.930000E+01 

1 • 700000E-0 1 

9 

1 5 

3.328300E-0 1 

1 . 6661 60E+05 

1.11 OOOOE-O 1 

15 

9 

3.4371 60E-0 1 

1 . 491 890E+05 

1 . 420000E— 0 1 

3 

1 

3.601 300E-0 1 

1 . 688930E+05 

1 • 340000E-0 1 

3 

9 

3 • 842700E-0 1 

1 . 706370E+05 

1 . 530000E-0 1 

9 

9 

3.919010E-01 

1 . 646 1 20E+05 

2. 31 0000E-01 

3 

3 

3.995000E-01 

1 . 49 1 890E+05 

6. 300000E— 0 1 

3 

5 

4.026080E-0 1 

1 . 86572 0E + 05 

2.800000E-01 

7 

9 

4.040900E-0 1 

1 . 865930E+05 

8 . 300000E-0 1 

21 

27 

4 . I 37400E-0 1 

2.056770E+05 

I . 1 70000E-0 I 

1 5 

5 

4.176160E-01 

1 . 870920E+ 05 

8 . OOOOOOE-O 1 

5 

7 

4.227750E-0 1 

1 . 7421 30E+05 

1 .71 OOOOE-O 1 

5 

3 

4 .239400E-0 1 

1 . 874720E+05 

8. 1 OOOOOE-O 1 

1 5 

2 1 

4.434600E-01 

1 . 888840E+05 

9. 1 OOOOOE— 0 1 

9 

1 5 

4.447030E— 0 1 

1 . 646 1 20E+05 

6.420000E— 01 

3 

5 

4.530400E— 01 

1 . 89336 0E+ 05 

6.700000E-01 

7 

9 

4.552540E-0 1 

1 . 89336 0E+05 

3.050000E-01 

7 

9 

4 . 623200E— 0 1 

1 .4901 30E+05 

3.370000E— 01 

9 

9 

4 .677930E-0 1 

1.90121 0E+05 

9. OOOOOOE-O 1 

3 

5 

4 . 793500E-0 1 

1 .6661 60E+05 

1 • 230000E-0 1 

1 5 

1 5 

5 • 000800E— 0 1 

1 . 688930E+05 

8.400000E— 01 

3 

9 

5.004000E-0 1 

1 . 666 1 60E+05 

6.390000E— 01 

15 

21 

5.006900E-0 1 

2 , 05677 0E+ 05 

2.91 0000E-01 

15 

1 5 

5.028800E-0 1 

1 .490 1 30E+O5 

9.700000E— 02 

9 

3 

5. 1 04450E-0 1 

1 • 78274 0E+05 

2. 220000E-0 1 

1 

3 

5. 1 75900E-0 1 

2.256440E+05 

5.600000E-01 

1 5 

25 

5. 1 77800E-0 1 

2.23731 0E+05 

5.700000E-01 

25 

35 

5 . 335800E-0 1 

2 . 25644 0E + 05 

1 . 780000E-0 1 

15 

1 5 

5.478800E-0 1 

1 .706370E+05 

1 • 800000E-0 1 

9 

9 

5.537400E-0 1 

2.056770E+05 

4 . 32 OOOOE-O 1 

15 

25 

5.679400E-0 1 

1 . 490 1 30E + 05 

4.520000E— 01 

9 

1 5 

5.938500E-0 1 

1 . 706370E + 05 

4.980000E-0! 

9 

I 5 

6.1681 OOE-O 1 

1 • 865930E + 05 

1 • 480000E-0 1 

21 

1 5 

6.242520E-0 1 

1 . 893360E + 05 

1 .420000E— 01 

7 

5 

6.345800E-0 1 

1 • 874720E + 05 

1.11 OOOOE-O 1 

1 5 

9 

6.482070E-0 1 

1 .491 890E + 05 

2.300000E-01 

3 

3 

6.61 0580E-0 1 

1 . 7421 30E + 05 

5. 400000E-0 1 

5 

7 

6.629800E-0 1 

1 . 870920E + 05 

1 . 1 20000E-0 1 

5 

3 


N2 1 4 l.OOOOOOE+OO 

0 6 3 

1. O.OOOOOOE+OO 

1. 7.223639E+04 

2. 8.578989E+04 

2. 9.969445E+04 

2. 1.282647E+05 
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1. 1.503260E+05 

3.395007E+03 2.079929E+01 2.891988E+00 2.690556E-02 

2*1011 80E+0 3 1 * 998637E + 0 1 2.071782E + 00 1.870440E-02 

2.495037E+03 2.081944E+01 2.356754E+00 2.647392E-02 

2.434464E+03 1.840369E+01 2.355316E+00 3.165360E-02 

2.928102E+03 2.457470E+01 2.627105E+00 2.834436E-02 

1.081460E+03 6.935016E+00 1.660375E+O0 6.906240E-03 

N2 FIRST POSITIVE BAND SYSTEM 
14 14 32 -1 3.000000E-03 

8.578989E+04 2.495037E+03 2.081944E+01 0.000000E+00 2.356754E+00 

7.223639E+04 2.101180E+03 1.998637E+01 -3 . 597000E-02 2.071 782E+00 


3382+00 

4064+00 

1 975+00 

5014-01 

7190-02 

5871-03 

2616-04 

5698-06 

4722-08 

6508-1 1 

2780-14 

81 14-15 

2515-15 

5962-16 

3248+00 

2310-02 

2120+00 

2987+00 

1318+00 

2729-0 1 

2925-02 

1613-03 

4174-05 

3963-07 

6055-10 

3322-13 

5035-15 

1662-1 5 

1 900+00 

1 032+00 

1 1 32+00 

3868-01 

2738+00 

2106+00 

6148-01 

8462-02 

5675-03 

1719-04 

1 848-06 

3035-09 

2138-12 

4334-15 

8857-01 

1 782+00 

1205-02 

1623+00 

1806-02 

1 808+00 

2605+00 

1065+00 

1 857-0 1 

1 495-02 

524 1 -04 

6322-06 

1101 -08 

1 051-1 1 

3649-01 

1 450+00 

7724-01 

3227-01 

1 139+00 

4780-01 

8305-0 1 

2706+00 

1 561+00 

3420-01 

3274-02 

1318-03 

1769-05 

3236-08 

1399-01 

8647-01 

1 275+00 

9050-02 

8823-01 

4262-01 

1 040+00 

1916-02 

2438+00 

2029+00 

5569-01 

6298-02 

2890-03 

4290-05 

51 46-02 

4367-01 

1 127+00 

6910-0 1 

5227-02 

1 057+00 

3171-02 

1290+00 

2766-04 

191 8+00 

2402+00 

8265-01 

1 099-0 1 

571 7-03 

1351-02 

2000-01 

7496-01 

1 009+00 

1798-01 

3829-01 

8078-0 1 

6748-02 

1 1 59+00 

1 694-0 1 

1 299+00 

2530+00 

1 1 40 + 00 

1 777-0 1 

6595-03 

8626-02 

4247-01 

9079-01 

6361-01 

1 230-05 

6967-01 

3946-01 

3633-01 

7862-01 

5219-01 

7206-01 

2688+00 

1 480+00 

2348-03 

3589-02 

2182-01 

5406-01 

8349-01 

2454-0 1 

1289-01 

7733-01 

8360-02 

5773-01 

3764-01 

3779-01 

2857-01 

2675+00 

8403-04 

1463-02 

1 053-01 

3917-0 1 

751 1-01 

5820-01 

2460-02 

3833-01 

6077-01 

3194-03 

8370-01 

9012-02 

1 1 09-0 1 

4687-02 

3038-04 

5907-03 

4884-02 

2187-0 1 

5510-01 

7056-01 

281 0-0 1 

2738-02 

5721-01 

3305-01 

1 307-01 

7897-01 

2617-04 

1 1 56+00 

1 1 14-04 

2380-03 

2212-02 

1151-01 

3567-01 

6338-01 

5261-01 

6400-02 

1 852-01 

5945-01 

9576-02 

3527-01 

5846-0 1 

9342-02 

4147-05 

9627-04 

9884-03 

5824-02 

2128-01 

4795-01 

6039-01 

291 2-01 

8677-04 

3729-01 

4603-0 1 

6687-04 

5423-01 

3235-01 


4.800000E-01 2.000000E+00 

N2 SECOND POSITIVE BAND SYSTEM 
7753-1 6.000000E— 02 

1.282647E+05 2.928102E+03 2.457470E+01 -3 . 093420E+00 2.627105E+00 

8 • 578989E+04 2.495037E+03 2.081944E+01 O.OOOOOOE+OO 2.356754E+00 
4493+00 3899+00 1349+00 2363-01 2190-02 0000+00 0000+00 

3287+00 1868-01 3223+00 2515+00 6957-01 0000+00 0000+00 

1469+00 2038+00 3299-01 1630+00 3034+00 0000+00 0000+00 

5226-01 2003+00 5957-01 1181+00 4752-01 0000+00 0000+00 

1634-01 1124+00 1614+00 1841-02 1570+00 0000+00 0000+00 

4726-02 4839-01 1427-01 8891-01 1416-01 0000+00 0000+00 

1300-02 1790-01 8303-01 1345+00 2939-01 0000+00 0000+00 

2.600000E-01 5.500000E— 01 
N2 BIRGE-HOPFIELD BAND SYSTEM NO. 2 
14 14 6 1 1 4.300000E— 01 

1.503260E+05 1.081460E+03 6.935016E+00 O.OOOOOOE+OO 1.660375E+00 
O.OOOOOOE+OO 3 . 395007E + 03 2.079929E +01 -2 . 876600E-02 2.891 988E+00 
2030-06 2673-05 1800-04 8201-04 8162-03 8162-03 1975-02 

4176-02 7869-02 1343-01 2100-01 3040-01 4107-01 5211-01 

2878-05 3336-04 1961-03 7784-03 5708-02 5708-02 1168-01 
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2062-01 

3197-01 

4407-01 

5445-0 1 

6053-01 

6050-01 

5404-0 1 

21 1 3-04 

2141-03 

1091-02 

3720-02 

1942-01 

1942-01 

3272—0 l 

4651-01 

5638-01 

5832-01 

5085-01 

3613-01 

1919-01 

5814-02 

1063-03 

9324-03 

4070-02 

1 173-01 

4130-01 

4130-01 

5475-01 

5846-01 

4940-0 1 

31 1 0-01 

1 209-0 1 

1 074-02 

1929-02 

1 196-01 

4098-03 

3076-02 

1 132-01 

2695-01 

5878-01 

5878-01 

5620-0 1 

3822-01 

1533-01 

1 250-02 

2797-02 

1551-01 

2806-01 

3106-01 

1284-02 

8129-02 

2468-01 

4702-01 

5468-01 

5468-01 

3050-01 

6758-02 

4457-03 

1221-01 

2752-01 

3141-01 

2129-01 

691 0-02 

3396-02 

1 779-0 1 

4332-01 

6286-01 

2809-01 

2809-01 

3394-02 

3245-02 

2080-01 

3262-01 

2576-01 

91 03-02 

6593-04 

5677-02 

7757-02 

3287-0 1 

6166-01 

6290-0 1 

3065-02 

3065-02 

5377-02 

2616-01 

3340-01 

1858-01 

21 14-02 

2807-02 

1 620-01 

2439-01 

1657-01 

5180-0 1 

7053-01 

4333-0 1 

4927-02 

4927-02 

2847-0 1 

3350-01 

1376-01 

6984-04 

9006-02 

2344-01 

2239-01 

8604-02 

2782-01 

6982-01 

6261-01 

1557-01 

2758-01 

2758-01 

3513-0 1 

1227-01 

1 776-03 

1439-01 

2641-0 1 

1694-01 

2066-02 

2643-02 

4466-01 

8004-01 

3936-01 

2292-04 

3852-0 1 

3852-01 

1459-01 

2976-03 

1745-01 

2734-01 

1223-01 

3466-04 

8758-02 

2074-01 

6484-01 

7664-01 

1304-01 

9250-02 

2196-01 

2196-01 

5961-06 

1773-01 

2839-01 

9869-02 

4287-03 

1434-01 

2229-01 

1 069-0 1 

8563-01 

5880-0 1 

1031-03 

3278-0 1 

1637-02 

1637-02 

1441-01 

3046-01 

1 026-01 

1 053-02 

1801-01 

2177-01 

5835-02 

9240-03 

1 033+00 

3266-01 

9064-02 

4681 -0 1 

6908-02 

6908-02 

3231-0 1 

1414-01 

7642-03 

1 971-01 

2150-01 

3233-02 

3700-02 

1801-01 


9.000000E-02 1.600000E-01 

N2 MOLECULE FREE-FREE RADIATION 
0 51 5.045000E-02 1.000000E-01 3.593699E-02 

N2+ 1 3 2.027000E+14 

0 4 2 

1 . 0 . OOOOOOE+OO 

2. 1.295200E+04 

1. 3.663400E+04 

1. 9.297800E+04 

3.175700E+03 2.321650E+01 2.779800E+00 2.877600E-02 

2.736600E+03 2.145250E+01 2.503500E+00 0. OOOOOOE+OO 

3.481670E+03 3.336580E+01 2.997000E+00 2.806000E-02 

2.949540E+03 O.OOOOOOE+OO 2.374000E+00 7.194000E-02 

N2+ MEINEL BAND SYSTEM 
77211 1.200000E-02 

1.295200E+04 2.736600E+03 2.145250E+01 O.OOOOOOE+OO 2.503500E+00 

O.OOOOOOE+OO 3.1 75700E+03 2.321650E+01 -5 . 755200E-02 2.779800E+00 

4751+00 3255+00 1360+00 4526-01 1329-01 3624-02 0000+00 

3798+00 3115-01 2245+00 1990+00 1032+00 4145-01 0000+00 

1226+00 3358+00 2137-01 7974-01 1745+00 1395+00 0000+0 

2055-01 2368+00 1851+00 1048+00 7195-02 1078+00 0000+00 

1914-02 6236-01 2946+00 6274-01 1553+00 6918-02 0000+00 

9822-04 7946-02 1171+00 2929+00 5808-02 1509+00 0006+00 

2573-05 5163-03 1972-01 1738+00 2479+00 5609-02 0000+00 

5.000000E— 01 5. OOOOOOE+OO 

N2+ FIRST NEGATIVE BAND SYSTEM 
14 14 3 1 -1 3.500000E-02 

3.663400E+04 3.481670E+03 3.336580E+01 -7 . 733550E-0 1 2.997000E+00 

O.OOOOOOE+OO 3. 175700E+03 2.321650E+01 -5. 755200E-02 2.779800E+00 

6509+00 3014+00 4537-01 2248-03 1452-04 4634-06 9482-08 

6438-09 4200-12 1275-11 3949-13 2877-15 5259-15 7421-19 

2588+00 2226+00 4060+00 1056+00 6935-02 3986-04 3088-05 
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341 7-07 

5755-08 

3877-1 1 

1028-10 

7000-12 

2491-14 

9910-15 

7016-01 

2860+00 

5065-01 

4137+00 

1660+00 

1340-01 

5729-04 

1 1 32-04 

4951-07 

2659-07 

1710-09 

3877-10 

5560-1 1 

1 009-12 

1600-01 

1 324+00 

2290+00 

2100-02 

3792+00 

2205+00 

2069-01 

4930-04 

3006-04 

1 742-07 

8356-07 

1676-08 

7849-10 

2755-10 

3297-02 

4273-0 1 

1654+00 

1557+00 

6726-02 

3310+00 

2673+00 

2789-01 

1807-04 

6420-04 

3295-07 

1969-06 

91 84-08 

5702-1 0 

6342-03 

1 1 40-0 1 

71 13-01 

I 706+00 

9290-01 

2925-01 

2830+00 

3068+00 

3419-0 1 

1067-05 

1 161-03 

6715-06 

3593-06 

351 9-07 

1 155-03 

2700-02 

2362-01 

9451-01 

1569+00 

4815-01 

5331-0 1 

2414+00 

3401+00' 

3890-01 

8094-04 

1823-03 

3522-05 

4993-06 

2000-04 

5861-03 

6691-02 

3801-01 

1 096+00 

1 333+00 

2044-01 

7236-01 

2083+00 

3686+00 

4153-01 

3872-03 

2509-03 

1 172-04 

3280-05 

1 1 85-03 

1695-02 

1261-01 

5236-01 

1 161+00 

1 065+00 

5939-02 

8469-0 1 

1837+00 

3933+00 

4170-01 

1 087-02 

301 0-03 

5034-06 

2239-04 

3934-03 

3668-02 

2003-01 

6488-01 

1 150+00 

8084-01 

4743-03 

9065-01 

1 672+00 

4149+00 

3928-01 

2361-02 

7027-07 

3931-05 

8445-04 

9639-03 

6598-02 

2826-01 

7442-01 

1084+00 

5865-0 1 

5552-03 

9141-01 

1579+00 

4338+00 

3432-01 

8419-08 

6286-06 

1675-04 

2324-03 

1942-02 

1 043-01 

3652-01 

8051-01 

9821-01 

4075-01 

3604-02 

8818-01 

1554+00 

4496+00 

7471-09 

8755-07 

3030-05 

5156-04 

5207-03 

3406-02 

1496-01 

441 1-01 

8319-01 

8628-01 

2710-01 

7896-02 

8204-01 

1 593+00 

2725-10 

9533-08 

4831-06 

1 044-04 

1281-03 

l 007-02 

5378-02 

1991-01 

5053-01 

8288-01 

7393-01 

1 720-01 

1 238-01 

7384-0 1 


2.500000E-0 1 6.500000E-01 

N2 + MOLECULAR ION FREE-FREE RADIATION 
0 51 5.030000E-02 l.OOOOOOE-Ol 1.000000E+00 


NO 1 4 2 « 694000 E+ 1 6 

0 3 2 

2. 8.640000E+01 

1 • 6.325785E+04 

2. 6.608747E+04 

2.739270E+03 2.01 0000E+01 

3.41 1830E+03 2.083380E+01 

1 .493020E+03 1 .09392 0E + 01 

NO GAMMA BAND SYSTEM 
7 7 2 1-1 2.500000E-02 

6.325785E+04 3.41 1830E+03 

8.640000E+0 1 2. 739270E+03 

1653+00 3291+00 2906+00 

2636+00 1050+00 1460-01 


2377+00 
1606+00 
9080-01 
4560-01 
21 10-01 
1 .900000E-0 1 
NO BETA BAND 
7 7 3 1 1 

6.608747E+04 
8.640000E+0 1 
2264-04 
31 1 0-03 
2046-02 
8569-02 


8000-03 
7210-01 
1 349+00 
1340+00 
9870-0 1 


1 545+00 
7530-01 
6000-03 
3420-01 
8840-01 


3 . 50000 OE — O 1 
SYSTEM 

1 .000000E-02 
1 .493020E+03 
2. 739270E+03 


1851-03 
2089-02 
1093-0 1 
3487-01 


7921-03 

7278-02 

2975-01 

6959-01 


2. 4 5258 OE +00 
2.870690E+00 
1 . 6208 1 OE+OO 


2.083380E+01 
2.0 1 OOOOE+O 1 
1504+00 
1904+00 
5000-0 1 
3090-0 1 
2950-0 1 
0000+00 
0000+00 


1 . 093920E+0 1 
2.010000E+01 
2364-02 
1 749-01 
5445-0 1 
8817-01 


2.561 060E-02 
2 . 359630 E — 02 
2.1941 70E-02 


-4.028640E-0 1 2.870690E+00 

-1 .726560E-03 2.452580E+00 


5080-01 

2370+00 

3660-01 

1526+00 

0000+00 

0000+00 

0000+00 


1 160-01 
1588+00 
1 077+00 
0000+00 
oooo+oo 
oooo+oo 
0000+00 


3200-02 

3520-01 

0000+00 

0000+00 

oooo+oo 

0000+00 

0000+00 


1 .391320E-01 
-1 • 72 6 560 E— 03 
5530-02 1080 

3249-01 4959 

7450-0 1 7975 

7569-01 4184 


1 . 62081 OE+OO 
2.452580E+00 
-01 1834-01 

— 0 1 6444-01 

-01 6747-01 

-01 1072-01 
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2560-01 

7461-01 

9835-01 

7012-01 

2192-01 

3716-04 

1382-01 

5799-01 

1 1 03+00 

7784-01 

1483-01 

3280-02 

3385-01 

4922-0 1 

1 033+00 

1097+00 

2228-01 

5985-02 

4703-01 

4824-01 

1336-01 


2.100000E-01 5 • OOOOOOE— 0 1 

NITRIC OXIDE FREE-FREE RADIATION 
0 5! 5.050000E-02 1.000000E-01 3.065373E-02 

NO INFRARED VIBRATION-ROTATION BANDS 
2.450000E+00 1.90385OE+03 1.397000E+01 1.704600E+00 1 »780000E-02 

0 30 300 6.000000E+03 1.666666E+01 


2942-05 

5969-05 

9081-05 

1 230-04 

1565-04 

1915-04 

2280-04 

2657-04 

3041-04 

3425-04 

I 138-06 

3388-06 

6636-06 

1 092-05 

1628-05 

2275-05 

3028-05 

3884-05 

4826-05 

5840-05 

6820-09 

2834-08 

8956-08 

2510-07 

6278-07 

1 392-06 

2743-06 

4844-06 


7771-06 1147-05 

003 8.692000E+18 

OXYGEN ATOM FREE-FREE CONTINUUM 
0 51 5.020000E-02 5.000000E-02 1 .268524E-02 

OXYGEN ATOM FREE-BOUND CONTINUUM (HENRY FORMULA) 
0 3 14 5.010000E-02 1 .000000E-04 

9. 1.12231 1E+02 

5. 2.283140E+04 

1. 4 . 862253E+04 

9. 2.000000E+00 1.12231 1E+02 2.162665E+05 

9. 2 . 00000 OE+OO 1. 12231 1E+02 1.966252E+05 

5. 2.000000E+00 2. 283140E+04 2.162665E+05 

5. 2 . 00 000 OE+OO 2.283140E+04 1 .966252E+05 

1. 2.000000E+00 4 • 862253E+04 2.162665E+05 

9. 2.000000E+00 1. 12231 1E+02 1.580396E+05 

9. 3.000000E+00 1 .404738E+05 1 .580330E+05 

16. 4 . 00 000 OE+OO 1.481559E+05 1.580330E+05 

25. 5 . 00 000 OE+OO 1.517117E+05 1.580330E+05 

36. 6 . 00 OOO OE+OO 1 .536432E+05 1.580330E+05 

49. 7 • 00 000 OE+OO 1 .548078E+05 1 .580330E+05 

64. 8.000000E+00 1 .555637E+05 1.580330E+05 

81. 9 . OOOOO OE+OO 1 .560820E+05 1 .580330E+05 

100. 1 . OOOOOOE+O 1 1 .564527E+05 1.580330E+05 


1 .500 

4.311 

2.260 

1 .500 

4.378 

3.850 

1 .500 

4.800 

1 .950 

1 .500 

6.829 

4.640 

1.500 

5*124 

7.650 

1 .000 

2.661 

2.940 

1 .000 

1 .000 

0.000 

1 .000 

1 .000 

0.000 

1 .000 

1 .000 

0.000 

1 .000 

1 .000 

0.000 

1 .000 

1 .000 

0.000 

1 .000 

1 .000 

0.000 

1 .000 

1 .000 

0.000 

1 .000 

1 .000 

0.000 


OXYGEN ATOM LINE RADIATION 
0 25 0 1.580330E+05 1.600000E+01 

8.113700E-02 7.800000E+01 7.700000E-03 9 9 

8.784500E— 02 7.800000E+01 3.700000E-02 9 9 

9.3601 10E-02 1.586800E+04 1.500000E-02 5 7 

9 . 894600E— 02 7.800000E+01 5.600000E-02 9 15 

9.994940E— 02 I.586800E+04 3.500000E-02 5 3 


58 



APPENDIX B - Continued 


1 .026600E-0 1 7 « 800000E+0 1 1.000000E-02 9 15 

1.152160E-01 1.586800E+04 9.000000E-02 5 5 

1.217640E-01 3.379200E+04 1.300000E-01 1 3 

1.303500E-01 7.800000E+01 3.100000E-02 9 3 

6.157300E-01 8.662900E+04 6.640000E-02 15 25 

6.653780E-01 1.159180E+05 1.330000E-01 3 1 

7.156800E-01 1 .026620E+05 3.630000E-01 5 5 

7.477300E-01 1.139160E+05 5.700000E-01 9 15 

7.773400E-01 7.376800E+04 9.220000E-01 5 15 

7.886310E-01 1.159180E+05 5.750000E-01 3 5 

7.949300E-01 1 • 0 1 1 430E+05 4.950000E-01 15 2l 

7.989900E-01 8.863100E+04 4.600000E-01 9 15 

8.226800E-01 1.011430E+05 3.270000E-01 15 15 

8.446500E-01 7.679500E+04 8.980000E-01 3 9 

8.508630E-01 1.159180E+05 B.140000E-01 3 3 

8.820450E-01 1.026620E+05 4.260000E-01 5 7 

9.263900E-01 8.662900E+04 9.000000E-01 15 25 

1.128700E+90 8.863100E+04 7.500000E-01 9 15 

1.129900E+00 8.662900E+04 1.730000E-01 15 5 

1 • 3 1 6400E+00 8.863100E+04 1.630000E-01 9 3 

0+031 .621 OOOE+1 6 
0+ FREE-FREE CONTINUUM 

0 51 5.010000E-02 5.050000E-02 l.OOOOOOE+OO 

0+ DUMMY FREE-BOUND (PARTITION FUNCTION DATA) 

0 3 1 5.020000E-02 1 .000000E-04 

4. O.OOOOOOE+OO 

lO. 3 . 85840 1 E+04 
6. 5.822450E+04 

1. l.OOOOOOE+OO l.OOOOOOE+OO O.OOOOOOE+OO 

1.000 1.000 0.000 

OXYGEN ION LINE RADIATION 
0 50 1 4.075310E+05 1.600000E+01 

5.181 3 00 E— 02 4 • 04670 0E+ 04 4.500000E-02 6 6 

5. 1 55500E— 02 4.046700E+04 9.700000E-01 6 10 

5.393700E— 02 O.OOOOOOE+OO l.lOOOOOE-Ol 4 12 

5.550800E— 02 2.681700E+04 6.700000E-02 10 10 

6.005850E-02 4.046700E+04 3.900000E-02 6 10 

6.441480E-02 4.046700E+04 1.500000E-01 6 2 

6.732200E-02 4.046700E+04 6.300000E-02 6 6 

7.185400E— 92 2.681700E+04 2.500000E-01 10 10 

7.966610E— 02 4.046700E+04 7.000000E-02 6 10 

8.338900E-02 O.OOOOOOE+OO 4.300000E-01 4 12 

2.738000E-01 2.039420E+05 1.200000E-01 2 6 

3.133900E— 01 2.068950E+05 1.350000E-01 20 12 

3 . 2 72400E— O 1 2.287370E + 05 1.370000E-01 14 10 

3.292400E-01 2.084310E+05 1.380000E-01 12 12 

3.385900E— O 1 2.039420E+O5 9.600000E-01 2 6 

3 • 408400E— 9 1 2.299560E+95 1.400000E-01 10 10 

3.469100E-O1 2.11 6360E+05 1.350000E-01 10 6 

3.733900E-01 2.325130E+95 1.450000E-01 6 10 

3.735900E— O 1 1.854020E+05 1.230000E-01 12 4 

3.753700E— 01 2.121620E+05 1.680000E-01 4 12 

3.809300E-01 2.142100E+05 1.420000E-01 6 6 

3.867200E-O1 2.068950E+05 1.300000E-01 20 20 

3.91 4400E —01 2 . 06972 0E + 05 1.930000E-01 10 6 

3.966900E-01 1.890080E+95 3.430000E-01 6 6 

4 . 074800E— O 1 2.068950E+05 6.900000E-01 20 28 

4.093700E— 01 2.314290E+05 8.400000E-01 28 36 
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4.111 400E-0 1 

2.0843.1 0E + 05 

6. 300000E— 0 1 

12 

20 

4. 1 51 700E-0 1 

2.08431 0E+05 

2.61 OOOOE-O 1 

12 

1 2 

4.187900E-01 

2.287370E+05 

8.500000E— 01 

14 

18 

4.272300E-0 1 

2.26851 0E+05 

8.900000E— 01 

2 

6 

4.278000E-0 1 

2 . 327460E+05 

8. 1 OOOOOE— 0 1 

20 

28 

4.293800E-0 1 

2.3251 1 0E+05 

9, 100000E— 01 

12 

20 

4.31 4300E-0 1 

1 .854020E+05 

2.97000OE— 01 

12 

12 

4.314700E-01 

2 . 32e890E+05 

8.400000E— 01 

14 

18 

4.3257O0E-01 

2.3251 1 0E+05 

1 .700000E— 01 

6 

2 

4.349700E-0 1 

2.069720E+05 

2 . 960000E— 0 1 

10 

1 0 

4.385300E-0 1 

2.11 6360E+05 

1 . 240000E— 0 1 

1 0 

1 0 

4 .4 1 81 OOE-0 1 

1 . 890080E+05 

5.500000E— 01 

6 

1 0 

4.4461 OOE-O 1 

2. 287370E+05 

1 . 760000E— 0 1 

14 

1 4 

4.467200E-01 

2 . 453960E+05 

8.300000E— 01 

6 

18 

4.593200E-0 1 

2.069720E+05 

4.900000E-01 

10 

1 4 

4.651500E-01 

1 . 854020E+05 

5 • 500000E-0 1 

12 

20 

4.698000E-01 

2.3251 1 0E+05 

3. 470000E-0 1 

6 

6 

4.700800E-01 

2 • 299550E+05 

4. 1 OOOOOE— 01 

10 

1 4 

4.703900E-0 1 

2.11 6360E+05 

6.400000E-01 

10 

1 4 

4.8681 OOE-0 1 

2.3251 1 0E+05 

2.590000E-01 

6 

1 0 

4.91 3000E-0 1 

2. 121 620E+05 

7.300000E-01 

4 

12 

4.943200E— 01 

2. 1421 00E+05 

6.500000E-01 

6 

1 0 

5.191 1 OOE-0 1 

2. 1421 00E+05 

2. 1 30000E-0 1 

6 

6 

6.897500E— 0 1 

2.314290E+05 

1 . 700000E-0 1 

28 

20 


02 1 2 3.536000E+13 

0 5 1 

3. O.OOOOOOE+OO 

2. 1*1 3924 7E + 04 

1. 1 • 8985 1 3E+04 

3. 5.193451E+04 

3. 7.165469E+04 

2 • 2 73805E+0 3 1.737049E + 01 2.079999E + 00 2.271847E-02 

2*171 1 32E+03 1.726546E+01 2.052288E+00 2.460328E-02 

2*061 334E+03 2.0071 10E+01 2.014894E+00 2.614279E-02 

1*1 7S368E+0 3 3.237274E + 01 1.510728E+00 0.000000E+00 

1.007670E+03 1.151319E+01 1 . 1 78368E+00 1.582667E-02 

02 SCHUMANN-RUNGE BAND SYSTEM 
14 14 5 1 1 4.000000E— 02 

7.165469E+04 1.007670E+03 1.151319E+01 -5.399774E-0 1 1.178368E+00 

O.OOOOOOE+OO 2 . 273805E+03 1.737049E+01 7.855786E-02 2.079999E+00 


1418-08 

1 338-07 

6657-07 

2325-06 

6399-06 

1 477-05 

2972-05 

5352-05 

8786-05 

1333-04 

1 888-04 

2516-04 

3172-04 

3799-04 

4680-07 

4065-06 

1868-05 

6046-05 

1 546-04 

3328-04 

6264-04 

1 059-03 

1637-03 

2347-03 

3152-03 

3999-03 

4817-03 

5533-03 

7335-06 

5809-05 

2442-04 

7247-04 

1 705-03 

3387-03 

5903-03 

9269-03 

1 336-02 

1 792-02 

2261-02 

2703-02 

3082-02 

3365-02 

7269-05 

5189-04 

1 971-03 

5299-03 

1 133-02 

2049-02 

3262-02 

4692-02 

6214-02 

7686-02 

8973-02 

9968-02 

1 060-01 

1 085-0 1 

51 1 7-04 

3247-03 

1 098-02 

2632-02 

5025-02 

8131-02 

1 160-01 

1 497-01 

1 783-0 1 

1 988-01 

2097-0 1 

2111-01 

2043-0 1 

1909-01 

2725-03 

151 1-02 

4464-02 

9348-02 

1 557-01 

2195-01 

2724-0 1 

3054-01 

3152-01 

3037-0 1 

2762-01 

2391-01 

1 984-01 

1587-01 

1 1 42-02 

5416-02 

1 364-01 

2423-0 1 

3402-01 

4010-01 

4 1 1 5-0 1 

3762-01 

3108-01 

2338-0 1 

1605-01 

1 000-01 

5574-02 

2688-02 

3865-02 

1527-01 

31 70-01 

4576-01 

51 18-01 

4669-01 

3551-01 

2245-01 

1 134-01 

4049-02 

6138-03 

6403-04 

1 146-02 

2790-02 

1077-01 

3419-01 

5571-01 

6090-01 

4862-01 

2844-01 

1098-01 
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1 650—02 

2457-03 

3576-02 

8194-02 

1 182-01 

1 356-01 

1 351-0 1 

2507-01 

6087-01 

7208-01 

5200-01 

2204-01 

3195-02 

7186-03 

8008-02 

1631-01 

2053-01 

1 990-0 1 

1609-01 

1 1 20-0 1 

6783-02 

4921-01 

8522-0 1 

6377-01 

2094-0 1 

4144-03 

6729-02 

2086-0 1 

2744-01 

2383-01 

1 505-0 1 

6688-02 

1620-02 

1 202-04 

6704-03 

821 7-01 

9095-0 1 

3108-01 

1 193-03 

1397-01 

3223-01 

3134-0 1 

1743-01 

4693-02 

2771-04 

2165-02 

6786-02 

1 049-0 1 

1 195-01 

1 1 74-00 

6854-0 1 

2393-02 

1685-01 

4101-01 

3129-01 

9194-02 

1 1 09-04 

5175-02 

1 389-0 1 

1812-0 1 

1 656-01 

1 1 79-01 

6697-02 

1442-00 

2905-0 1 

9461-02 

4830-0 1 

3394-01 

4514-02 

2764-02 

1697-01 

2421-01 

1930-01 

9601-02 

2400-02 

1 062-04 

1 101-02 


1 •700000E-0 1 3. 80000 0E- 01 

OXYGEN MOLECULE FREE-FREE CONTINUUM 
O 51 5«0 1 OOOOE— 02 5.000000E-02 2*540000E“02 
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TABLE I.- PARAMETERS FOR CALCULATION OF NEUTRAL BREMSST RAH LUNG 


obtained from eq. (12) for species H 2 , CO, C 2 , NO, and CN, using polarizability 
tables in ref. 37, p. 950. All other values of obtained from eq. (10) and 

calculations of refs. 12 and 13] 


Species 

Z^ used in program 

C 

0.02473 

O 

.01269 

H 

.07315 

H e 

.02256 

H 2 

.0107 

CO 

.0167 

°2 

.0174 

c 2 

.01535 

N 

.018 

NO 

.020 

N 2 

.0361 

CN 

.0167 



TABLE H.- NEUTRAL BREMSSTRAHLUNG CROSS SECTIONS (CARBON ATOM) 


Values in parentheses are taken from ref. 12. Values not in parentheses are 
computed from eqs. (8) and (10) for = 0.02473] 


|i m 

* 10 40 , cm^ 

T = 2000 K 

T = 5000 K 

. 

T = 10 000 K 

T = 15 000 K 

0.5 

5.97 

4.16 

5.52 

8.91 


(2.48) 

(4.45) 

(7.10) 

(9.67) 

1 

8.71 

11.04 

26.61 

47.14 


(5.29) 

(12.3) 

(26.1) 

(42.4) 

2 

18.54 

53.23 

139.0 

232.7 


(16.8) 

(51.9) 

(139) 

(254) 

5 

133.1 

463.6 

1064 

1673 


(129) 

(539) 

(1720) 

(3380) 
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TABLE HI.- PHOTOIONIZATION CROSS SECTIONS 


[The power-of-ten multiplier appears as a signed integer at the right of each entry 
(for example, 7.014 - 18 = 7.014 X 10“^®J 


X, jum 


0.0567 

.0700 

.0833 

.0967 

.1100 

.1101 

.1132 

.1162 

.1193 

.1223 

.1239 

.1240 

.1285 

.1330 

.1376 

.1421 

.1444 

.1445 

.3324 

.5204 

.7084 

.8963 

.9903 

.9904 

.1616 

.3327 

.5039 

.6750 

.7606 


T = 5000 K 

7.014-18 

9.364- 18 
1.139-17 
1.281-17 
1.333-17 
3.367-19 
3.408-19 
3.442-19 
3.471-19 
3.493-19 
3.501-19 
2.580-21 
2.609-21 
2.628-21 
2.638-21 
2.638-21 
2.634-21 
9.312-29 
1.135-27 
4.354-27 
1.098-26 
2.225-26 
3.001-26 
9.705-27 
1.566-26 

2.365- 26 
3.398-26 
4.695-26 
5.452-26 


(a) Carbon 


a FB , cm 2 


T = 8000 K 

T = 12 000 K 

6.932-18 

6.875-18 

9.208-18 

9.117-18 

1.117-17 

1.105-17 

1.256-17 

1.243-17 

1.311-17 

1.299-17 

1.355-18 

1.682-18 

1.372-18 

1.702-18 

1.386-18 

1.720-18 

1.397-18 

1.734-18 

1.406-18 

1.746-18 

1.410-18 

1.750-18 

5.239-20 

8.463-20 

5.297-20 

8.558-20 

5.336-20 

8.622-20 

5.356-20 

8.655-20 

5.357-20 

8.657-20 

5.349-20 

8.645-20 

1.402-23 

1.024-22 

1.709-22 

1.249-21 

6.555-22 

4.790-21 

1.653-21 

1.208-20 

3.349-21 

2.448-20 

4.517-21 

3.301-20 

2.476-21 

1.941-20 

3.995-21 

3.130-20 

6.033-21 

4.728-20 

8.669-21 

6.793-20 

1.198-20 

9.387-20 

1.391-20 

1.090-19 


TABLE HI.- PHOTOIONIZATION CROSS SECTIONS - Continued 


(b) Nitrogen 


X , jum 

ct fb , cm 2 

T = 5000 K 

T = 8000 K 

T = 12 000 K 

0.0838 

1.253-17 

1.246-17 

1.241-17 

.0844 

1.246-17 

1.240-17 

1.236-17 

.0850 

1.238-17 

1.234-17 

1.231-17 

.0852 

1.235-17 

1.231-17 

1.228-17 

.0853 

1.247-19 

1.853-18 

2.755-18 

.0860 

1.246-19 

1.853-18 

2.754-18 

.0866 

1.245-19 

1.852-18 

2.753-18 

.0873 

1.244-19 

1.850-18 

2.751-18 

.0879 

1.242-19 

1.848-18 

2.748-18 

.0883 

1.241-19 

1.846-18 

2.746-18 

.0884 

6.208-20 

1.004-18 

1.529-18 

.0902 

6.241-20 

1.009-18 

1.538-18 

.0920 

6.267-20 

1.014-18 

1.545-18 

.0937 

6.285-20 

1.017-18 

1.550-18 

.0955 

6.296-20 

1.019-18 

1.553-18 

.0964 

6.298-20 

1.020-18 

1.554-18 

.0965 

6.129-20 

9.278-19 

1.387-18 

.0977 

6.130-20 

9.282-19 

1.388-18 

.0989 

6.127-20 

9.282-19 

1.388-18 

.1002 

6.120-20 

i 

9.276-19 

1.387-18 

.1014 

6.111-20 

9.265-19 

1.386-18 

.1020 

6.105-20 

9.257-19 

1.385-18 

.1021 

2.349-21 

1.274-19 

2.318-19 

.1045 

2.357-21 

1.278-19 

2.325-19 

.1069 

2.360-21 

1.280-19 

2.328-19 

.1094 

2.359-21 

1.280-19 

2.327-19 

.1118 

2.354-21 

1.277-19 

2.322-19 

.1130 

2.349-21 

1.274-19 

2.317-19 

.1131 

1.750-31 

7.233-25 

8.885-24 

.2586 

2.091-30 

8.639-24 

1.061-22 

.4041 

7.976-30 

3.296-23 

4.049-22 

.5496 

2.006-29 

8.291-23 

1.018-21 

.6951 

4.058-29 

1.677-22 

2.060-21 



TABLE HI.- PHOTOIONIZATION CROSS SECTIONS - Concluded 


(c) Oxygen 


X, (im 

a FB , cm 2 

T = 5000 K 

T = 8000 K 

T = 12 000 K 

0.0710 

8.391-18 

8.494-18 

8.544-18 

.0725 

8.338-18 

8.436-18 

8.484-18 

.0732 

8.306-18 

8.402-18 

8.449-18 

.0733 

3.689-18 

4.012-18 

4.169-18 

.0736 

3.693-18 

4.015-18 

4.171-18 

.0738 

3.698-18 

4.018-18 

4.174-18 

.0740 

3.703-18 

4.021-18 

4.176-18 

.0743 

3.707-18 

4.023-18 

4.178-18 

.0744 

3.709-18 

4.025-18 

4.178-18 

.0745 

3.694-18 

3.866-18 

3.951-18 

.0763 

3.728-18 

3.887-18 

3.966-18 

.0782 

3.761-18 

3.904-18 

3.976-18 

.0800 

3.791-18 

3.918-18 

3.982-18 

.0819 

3.819-18 

3.929-18 

3.984-18 

.0828 

3.832-18 

3.932-18 

3.983-18 

.0829 

3.800-18 

3.620-18 

3.538-18 

.0835 

3.810-18 

3.629-18 

3.547-18 

.0842 

3.819-18 

3.639-18 

3.556-18 

.0848 

3.829-18 

3.647-18 

3.564-18 

.0855 

3.838-18 

3.656-18 

3.573-18 

.0858 

3.843-18 

3.660-18 

3.577-18 

.0859 

3.844-18 

3.654-18 

3.562-18 

.0871 

3.859-18 

3.669-18 

3.576-18 

.0882 

3.874-18 

3.683-18 

3.590-18 

.0894 

3.888-18 

3.696-18 

3.603-18 

.0905 

3.902-18 

3.709-18 

3.615-18 

.0911 

3.908-18 

3.715-18 

3.622-18 

.0912 

2.507-31 

4.451-25 

4.978-24 

.2530 

5.353-30 

9.505-24 

1.063-22 

.4148 

2.359-29 

4.189-23 

4.685-22 

.5767 

6.337-29 

1.125-22 

1.258-21 

.7385 

1.331-28 

2.363-22 

2.643-21 


) 


TABLE IV.- MOLECULAR BAND SYSTEMS USED IN PROGRAM SPRAD 


Species 

Band 

W ^ 

f g£ Ref. no. 

Franck- Condon 
factor 
(Ref. no.) 

H 2 

Lyman 

0.1103 

0.2 

1 

38 


Werner 

.1006 

.4 

1 

39 

02 

Schumann-Runge 

.2026 

.04 

26 

28 

n 2 + 

Meinel 

1.1091 

.012 

28 

28 


First negative 

.3914 

.035 

28 

28 

N 2 

First positive 

1.0510 

.003 

28 

28 


Second positive 

.3371 

.06 

28 

28 


Birge-Hopfield No. 2 

.0964 

.43 

28 

28 

NO 

Gamma 

.2263 

.025 

28 

28 


Beta 

.2199 

.01 

28 

28 

CO 

Angstrom 

.4511 

.003 

40 

26 


Hopfield-Birge (B-X) 

.1150 

.012 

41 

26 


Fourth positive 

.1544 

.15 

42 

41 

C 2 

Fox-Herzberg 

.2512 

.818 

43 

26 


Swan 

.5165 

.013 

44 

26 

o 

o 

+ 

First negative 

.2190 

.017 

40 

41 


Comet -tail 

.4910 

.0022 

41 

41 
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TABLE V.- THERMODYNAMIC STATE AND RADIATIVE CONTRIBUTORS 

IN SHOCKED C0 2 

Shock speed; 10.97 km/sec (36 000 ft/sec); 
initial pressure, 0.0133 N/cm2 (1 torr); 
initial temperature, 300 K; 
mixture temperature, 11 168 K; 
mixture pressure, 26.64 N/cm2 (2.63 atm)J 


Species 

e“ 

Concentration 
1.209E+ 17 

Contributing radiation 

CO 

2.241E+ 15 

(4+), Angstrom (B-A), Hopfield-Birge (B-X), free-free 

°2 

3.536E+ 13 

Schumann-Runge, free-free 

c 

4.329E+ 17 

Free-free, free-bound, lines 

0 

1.055E+ 18 

Free-free, free-bound, lines 

c + 

9.794E+ 16 

Free-free, lines 

o + 

1.621E+ 16 

Free-free, lines 

c ++ 

3.337E+ 15 

Free-free, lines 

C 2 

3.873E+ 13 

Swan, Fox-Herzberg (B-X), free-free 

co + 

1.192E+ 14 

Comet-tail, (1-), free-free 
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Figure 1.- Brief flow chart of spectral properties program SPRAD. 







Wave Number, l/cm 


Figure 3.- Emissivity of nitric oxide NO computed from smeared-line model of reference 41 

(N Mn = 3.215(10) 19 , T = 5000 K, L = 5 cm. 
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